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ABSTRACT 

A nonlinear  mathematical  model  has  been  formulated  of  a 
craft  having  a constant  dcadrisc  angle,  planing  in  regular  waves, 
using  a modified  low-aspect-ratio  or  strip  theory.  It  was 
assumed  that  the  wavelengths  would  be  large  in  comparison  to 
the  craft  length  and  that  the  wave  slopes  would  be  small.  The 
coefficients  in  the  equations  of  motion  were  determined  by  a 
combination  of  theoretical  and  empirical  relationships.  A 
simplified  version  for  the  case  of  a craft  or  model  being  towed 
at  constant  speed  was  programed  for  computations  on  a digital 
computer,  and  the  results  were  compared  with  existing  experi- 
mental data.  Comparison  of  computed  pitch  and  heave 
motions  and  phase  angles  with  corresponding  experimental  data 
was  remarkably  good.  Comparison  of  bow  and  center  of 
gravity  vertical  accelerations  was  fair  to  good. 

ADMINISTRATIVE  INFORMATION 

This  investigation  was  authorized  by  the  Naval  Sea  Systems  Command  with  initial 
funding  under  Task  Area  .SR;023-0 101  and  completion  under  Task  Area  ZF-43-42100], 

INTRODUCTION 

Computer  programs  for  estimating  the  motions  of  displacement  ships  in  waves  for  all 
headings  and  speeds  have  been  in  existence  for  some  time.  Comparable  computational 
schemes  for  planing  craft  do  not  exist  except  in  limited  and  restricted  cases.  A program  for 
planing  craft  would  be  quite  useful  to  the  small  craft  designer,  providing  a means  for 
systematically  exploring  the  effects  of  numerous  design  variations  on  performance  of  the 
craft  in  waves.  With  minor  modification,  the  program  could  also  be  used  to  examine  the 
merits  of  a hybrid  craft  design,  e.g.,  a combination  of  planing  craft  and  hydrofoil. 

Predicting  the  motions  of  a planing  craft  in  wave’s  is  by  no  means  a simple  problem. 
The  analytical  description  of  a high-speed  craft,  planing  in  waves,  involves  several  different 
types  of  flow  phenomena,  including  planing;  hydrodynamic  impact,  and,  to  a lesser  extent, 
surface  wave  generation  and  hydrostatics.  Also,  the  mathematics  tend  to  become  nonlinear 
rapidly  as  the  motion  increases  or,  like  the  real  craft,  can  in  some  instances  exhibit  large 
instabilities  such  as  porpoising. 

Development  of  a computer  program  that  would  take  into  account  all  of  the  previously 
described  factors  and  would  be  applicable  for  a wide  range  of  speed  and  wave  conditions 
requires  a careful  and  systematic  study  in  sevcial  stages  with  appropriate  verification  at  each 
stage.  To  lay  the  foundation  for  such  a general  program,  a simpler  problem  has  been 


formulated  in  this  report  with  potential  for  expansion  and  generalization  to  the  more 
complicated  case.  The  simpler  problem  is  that  of  a V-shaped  prismatic  body  with  hard  chines 
and  constant  deadrise  planing  at  high  speed  in  regular  head  waves. 

The  mathematical  formulation  is  analogous  to  low-aspect-ratio  wing  theory  with 
provisions  for  including  hydrodynamic  impact  loads,  essentially  a strip  theory.  Surface  wave 
generation  and  forces  associated  with  unsteady  circulatory  flow  are  neglected,  and  the  flow  is 
treated  as  quasi-steady.  The  mathematical  formulation  is  an  empirical  synthesis  of  several 
theoretically  derived  flows  describing  the  overall  craft  hydrodynamics.  Wave  input  is  restricted 
to  monochromatic  linear  deepwater  waves  with  moderate  wavelengths  and  low  wave  slopes. 

MATHEMATICAL  FORMULATION 

GENERAL 

Consider  a fixed  coordinate  system  (x,z)  (Figure  1)  with  x axis  in  the  undisturbed  free 
surface,  pointing  in  the  direction  of  craft  travel,  and  the  z axis,  pointing  downward.  If  the 
motions  of  the  craft  are  restricted  to  pitch  Q,  heave  and  surge  x^.^,  the  equation  of 
motions  can  be  written  as 

Mx^,q  = T^  - N sin  0 - D cos  0 

MZ(,q  = T^  - N cos  0 + D sin  0 + W 

10  =Nx^-Dxd+TXp  (1) 


where  M 
1 

N 

D 

W 

Tx 

T, 


is  mass  of  craft 

is  pitch  moment  of  inertia  of  craft 
is  hydrodynamic  normal  force 
is  friction  drag 
is  weight  of  craft 

is  thrust  component  in  x direction 
is  thrust  component  in  z direction 


X 


c 


distance  from  center  of  gravity  (CG)  to  center  of  pressure  for  normal  force 


Xj  is  distance  from  CG  to  center  of  action  for  friction  drag  force 
Xp  is  moment  arm  of  thrust  about  CG. 


Equation  (I)  is  exact;  however,  defining  the  hydrodynamic  forces  and  moment®  iii  waves 
can  be  extremely  difficult. 
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A high-speed  craft  moving  in  waves  may  transit  through  several  regimes  that  have 
different  hydrodynamic  flow  characteristics.  For  example,  as  the  craft  moves  away  from  the 
crest  of  wave,  the  flow  may  be  characterized  by  unsteady-state  planing  until  the  craft  collides 
with  the  oncoming  wave  crest  and  enters  another  regime  in  which  impact  forces  arc  important. 
After  the  impact,  the  craft  may  enter  still  another  regime  in  which  it  is  planing  but  in  which 
buoyancy  forces  are  rather  significant. 

The  most  promising  approach  to  a method  that  would  incorporate  all  three  types  of  flow 
conditions  into  a general  formulation  would  seem  to  be  a modified  strip  theory.  The 
mathematical  justification  for  this  approach  is  not  rigorous;  however,  there  is  sufficient 
precedent  to  expect  promising  results.  For  example,  impact  loads  on  landing  seaplanes  can 
be  estimated  reasonably  well  using  a strip  theory  incorporating  the  Wagner  two-dimensional 
(2-D),  expanding-wedge  theory,'  and  Chuang^  has  provided  a strip  method  for  determining  loads 
on  an  impacting  prismatic  form  that  agrees  extremely  well  with  experimental  results. 

More  recently,  Martin^  has  developed  a linear  strip  theory  for  estimating  motions  of  a 
planing  craft  at  high  speed,  which  shows  good  agreement  with  experimental  results.  A 
nonlinear  model  of  the  equations  of  motion  would  be  expected  to  provide,  in  addition  to  the 
motions,  reasonable  estimates  of  the  vertical  accelerations  which  are  an  important  consideration 
in  designing  a planing  craft. 

TWO-DIMENSIONAL  HYDRODYNAMIC  FORCE 

Implicit  with  any  strip  method  is  the  need  to  define  the  2-D  hydrodynamic  force  acting 
on  an  arbitrary  cross  section  of  the  body.  The  2-D  flow  problem  is  not  simple;  however,  it 
lends  itself  to  an  empirical  approach,  using  a combination  of  techniques  used  in  hydrodynamic 
impact  and  low-aspect-ratio  theories. 

The  typical  cross  section  of  a hard-chine,  V-shaped  prismatic  body  such  as  that  being 
considered  here  is  shown  in  Figure  2.  Figure  2 actually  illustrates  two  different  idealized- 
flow  conditions,  assumed  to  represent  the  crossflow  during  unsteady  planing,  depending  upon 
whether  the  flow  separates  from  the  chine  (Figure  2a)  or  not  (Figure  2b).  Nonwetted-chine 
flow  conditions  are  typical  of  the  sections  near  the  leading  edge  of  the  wetted  length  of  the 
craft.  Wetted-chine  flow  conditions  are  more  typical  of  sections  near  the  stern,  except 
possibly  in  the  most  extreme  motion  and  wave  conditions.  Some  sections  between  leading 
edge  and  stern  may  alternate  between  flow  conditions  as  the  wetted  length  changes  with  the 
motions. 


*A  complete  listing  of  references  is  given  on  page  33. 
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The  normal  hydrodynamic  force  per  unit  length  f,  acting  at  a section,  is  treated  as 
quasi-steady  and  is  assumed  to  contain  components  proportional  to  the  rate  of  change  of 
momentum  and  the  velocity  squared  (drag  term),  i.c. 

r=-(^(m,V)  + C„,,pbV=) 

where  V is  the  velocity  in  plane  of  the  cross  section  normal  to  the  baseline 
m^  is  the  added  mass  associated  with  the  section  form 

Cjj  ^ is  the  crossflow  drag  coefficient 

p is  the  density  of  the  fluid 
b is  the  half  beam. 

For  sections  near  the  leading  edge  of  the  wetted  length  with  nonwetted  chine,  the 
added  mass  Is  assumed  to  be  defined  in  the  same  manner  as  during  an  impact  which  for  a 
V-shaped  wedge  Is  given  by 


m3  = k3  7r/2pb^  (3) 

where  kg  is  an  added-mass  coefficient  that  may  also  include  a correction  for  water  pileup- 
kg  is  assumed  to  be  1.0  without  pilcup  correction. 

The  rate  of  change  of  momentum  of  the  fluid  at  a section  is  given  by 

§i  (niaV)  = mgV  + Vihg  - ~ (mg  V)  ^ (4) 

where  ^ is  the  body  coordinate  parallel  to  the  baseline;  see  Figure  1.  The  last  term  on  the 
right-hand  side  of  Fquation  (4)  takes  into  account  the  variation  of  the  section  added  mass 
along  the  hull.  This  contribution  can  be  visualized  by  considering  the  2-D  flow  plane  as  a 
substantive  surface  moving  past  the  body  with  velocity  U = -d^/dt  tangent  to  the  baseline. 

As  the  surfac ' moves  past  the  body,  the  section  geometry  in  the  moving  surface  may  change 

« 

with  a resultant  change  in  added  mass.  This  term  e.xists  even  in  steady-state  conditions  and 
IS  the  lift-producing  factor  in  low-aspect-ratio  theory. 

The  added  mass  of  a section  with  fully  wetted  chines  has  not  been  developed  to  the 
s.iine  ■-•xtenf  as  the  V wedge.  In  steady-state  planing  problems  such  as  those  of  Shuford,'^ 
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the  crossflow  is  treated  as  a Helmholtz-type  flow  in  which  the  Bobyleff  results  are  used  for 
estimating  drag  coefficients.  Helmholtz  flows  are  applicable  only  to  steady-state  conditions; 
so,  it  is  assumed  that  the  added  mass  for  the  fully  wetted  chine  flow  can  be  determined  from 
Equation  (3)  using  the  value  of  the  half-beam  at  the  chine.  In  using  the  Shuford  approach, 
it  is  assumed  that  the  crossflow  drag  i lefficient  for  a V-section  is  equal  to  the  drag  of  a flat 
plate  (Cp  g = 1.0)  corrected  by  the  Bobyleff  flow  coefficient  approximated  by  cos  (3,  i.e. 

Cq  c = 1.0  cos  |3  (5) 

The  Bobyleff  flow  coefficient  is  the  theoretical  ratio  of  the  pressure  on  a V-section  to  that 
experienced  by  a flat  plate  for  a Helmholtz-type  flow. 

The  same  approximation  is  used  for  estimating  the  drag  coefficient  for  nonwetted  chine 
sections,  using  the  instantaneous  value  of  the  half-beam  at  the  free  surface. 

An  additional  force  acting  on  the  body  is  the  buoyancy  force  fg.  This  force  is  assumed 
herein  to  act  in  the  vertical  direction  and  to  be  equal  to  the  equivalent  static  buoyancy  force 
multiplied  by  a correction  factor,  i.e. 


fg  = -apg(A)  (6) 

where  A is  the  cross-sectional  area  of  the  section,  and  a is  a correction  factor. 

The  full  amount  of  the  static  buoyancy  is  not  realized  because  at  planing  speeds  the  water 
separates  from  the  transom  and  chines,  reducing  the  pressure  at  these  locations  to  atmospheric 
or  less  than  the  equivalent  hydrostatic  pressure.  A greater  reduction  is  realized  in  the 
buoyancy  moment  because  of  the  corresponding  shift  in  the  center  of  pressure.  Shuford'* 
in  his  work  on  steady-state  planing  recommended  a factor  of  one-half  to  obtain  the  correct 
buoyancy  force.  In  the  following  computations,  the  buoyancy  force  was  corrected  by  a 
factor  of  one-half,  i.e.,  a = 1/2.  The  buoyancy  moment,  computed  as  the  static  buoyancy 
force  multiplied  by  its  corresponding  moment  arm,  was  corrected  by  an  additional  factor  of 
one-half  to  obtain  the  proper  mean-trim  angles. 

Equation  (2)  is  a synthesis  of  several  idealized  flow  conditions  combined  in  an  empirical 
manner.  In  all  of  these  flows,  it  is  assumed  that  the  net  relative  movement  of  the  fluid  past 
the  body  is  in  an  upward  direction.  This  condition  may  not  always  be  met  in  the  case  of 
unsteady  planing  in  waves.  Closer  scrutiny  will  be  required  to  determine  what  limitations 
will  be  imposed  upon  the  problem  as  formulated  and/or  what  modifications  will  be  required 
to  improve  the  formulation. 


5 


TOTAL  HYDRODYNAMIC  FORCE  AND  MOMENT 

The  total  normal  hydrodynamic  force  acting  on  the  body  is  obtained  by  integrating  the 
stripwise,  2*D,  hydrodynamic  force  given  by  Equations  (2)  and  (6)  over  the  wetted  length  C 
of  the  body.  A body  coordinate  system  {{.f)  with  its  origin  at  CG  and  the  | axis  pointing 
forward  parallel  to  the  baseline  of  the  body  is  defined  in  Figure  I to  facilitate  this  integration 
The  hydrodynamic  force  acting  in  the  vertical  or  z direction  of  the  fixed  integral  coordinate 
system  is  given  by 


N cos  6 


F,(t) 


■ [ cos  6 d{  + J d{ 


mg(t.t)V(tt)  + m,(«,t)V(tt) 


-U(t.ti  ^ [mj({.t)V({,t)l 
+ Cu  g(t.t)pba.t)V»{t,t))cosedt 


“I 


t apgAdI 


(7) 


where  the  integration  is  taken  over  the  Instantaneous  wetted  length,  Similarly  the  force 
acting  in  the  horizontal  or  x direction  is  given  by 


F^  ■ / fsinfld^ 


-U({,t)  ^ (mj(tt)V(|.t)l 
+ CD^^(«,t)pb(£.t)VJ(t.t))sln0dt 


(8) 


Wave  forces  are  obtained  by  neglecting  diffraction  and  assuming  that  the  wave  excitation 
is  caused  both  by  the  geometrical  properties  of  the  wave,  altering  the  wetted  length  and 
draft  of  the  craft,  and  by  the  vertical  component  of  the  wave  orbital  velocity  at  the  surface 
w^,  altering  the  normal  velocity  V.  The  horizontal  component  of  orbital  velocity  is  neglected. 
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since  it  is  assumed  small  in  comparison  with  the  forward  speed  The  velocities  U and  V 
may  then  be  written  as 

U » cos  0 - -w^)  sin  0 

V ■ x^.Q  sin  0 - + (Zj.Q-w^)  cos  0 (9) 

The  depth  of  submergence  h of  the  body  at  any  point  P({,f)  may  be  determined  by 

h ■ t sin  0 + f cos  0 - r (10) 

where  r is  the  instantaneous  vaiue  of  the  wave  elevation  directly  above  the  point. 

For  regular  head  waves  the  wave  elevation  for  a linear  deepwater  wave  is 

r ■ r^j  cos  k(x  + ct)  (II) 

where  rQ  Is  the  wave  amplitude 

k IS  the  wave  number 

c is  the  wave  celerity. 

At  point  P(t,?) 

X ■ XpQ  + t cos  0 + f sin  0 (1 2) 


(13) 


f nt.Dfdi-  f lb  cosOJdt 

A -'k 

“j^(m,(t,t)V(M)  + m,(tt)V(t.l) 

- U(tt)  ^ (m,(t.t)V({,t)) + Co_g(M)pb({,t)V2(t,t) 

+ apgA  cos  fij  (d{ 

EQUATIONS  OF  MOTION,  GENERAL 

Integrating  the  first  term  in  Equations  (7),  (8),  and  (13)  provides  hydrodynamic  forces 
and  moments  proportional  to  acceleration  of  the  motion.  These  can  be  combined  with  the 
inertial  terms  of  the  rigid  body  to  give  the  following  equation  of  motion 

(M  + M,  sin^  ^)XcG  ® cos  0)j!(,Q  - (Q^  sin  0)9 

■ T^ + F;, -Dcosfl  (14) 

(M,  sin  9 cosd)Xj.Q  +(M+M,  cos^  ®)*2cg 

- T^  + F;  + D sin  d + W 

-(Q,  sin  0))(^Q  -(Qj  cos  0)*Zj.(.  +(l  + lj)9 

■ Ffl  -Dxj  +TXp 


where  Mj(t)  •y 

f m,({,t)dt 

Q,(t) 

r mg{M){d« 

8 

'a(‘^  ■; 

f m,(t,t)«2d« 

fi 


F'j^  ■ F*  ■ ("(Mj  sin^  ^ 1 

Fy  ■ " {appropriate  acceleration  terms} 

F^  “ Fg  - (appropriate  acceleration  terms). 

A detailed  evaluation  of  the  integral  exprc.ssions  for  the  hydrodynamic  forces  and 
moments  is  provided  In  Appendix  A. 
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The  solution  to  Etiuation  (14)  is  cumbersome;  however,  it  can  be  accomplished  using 
standard  numerical  techniques.  Introducing  the  state  vector  (xj,  Xj,  Xj,  x^,  Xj,  Xj) 
where  x,  ■ 

’‘2  • *CC 

X3  -0 

^ “ ’‘CG 
’‘S  " *CG 

Xg  -0 

Equation  (14)  can  be  rewritten,  using  matrix  algebra,  as 

aT  ■ "g  (15) 

so  that 


T-A*'g  (16) 

where  A"’  is  inverse  of  the  inertial  matrix  A.  Equation  (16)  is  now  in  a form  that  lends 
Itself  to  Integration  by  using  a numerical  method  such  as  the  Runge>Kut>a>Merson  Integration 
routine. 

EQUATIONS  OF  MOTION,  SIMPLIFIED  FOR  CONSTANT  SPEED 

Assuming  that  the  perturbation  velocities  in  the  forward  direction  are  small  in  comparison 
to  the  speed  of  the  craft,  the  equations  of  motion  may  be  further  simplified  by  neglecting 
the  perturbations  and  setting  the  forward  velocity  equal  to  a constant,  i.e. 

- CONSTANT 

If  it  is  also  assumed  that  the  thrust  and  drag  forces  are  small  in  comparison  to  the  hvdrody- 
namic  forces  and  that  they  are  acting  through  the  center  of  gravity,  the  equations  of  mot'on 
may  be  written  as 


(M  + Mg  cos^  ~ ~ Fz 

-(Qj  cosfl)Zj.(,  +(l  + lg)0  = Fg 

These  equations  also  represent  the  ease  of  the  craft  (model)  being  towed  through  CG  at 
CONSTANT  speed.  Based  upon  the  previously  described  equations  of  motion,  a computer 
program  has  been  written  in  FORTRAN  language  to  compute  the  motions  of  a prismatic  body, 
planing  In  regular  head  waves  at  high  speed.  A listing  of  the  program  along  with  the 
appropriate  flow  chart  is  presented  in  Appendix  B.  The  listing  contains  reference  to  thrust 
and  drag  terms;  however,  they  have  no  significance,  except  to  provide  a starting  point  for 
possible  updating  of  the  program  to  Include  these  terms  in  the  future. 

COMPARISON  OF  COMPUTED  RESULTS  WITH  EXPERIMENTS 

Computations  of  pitch  and  heave  motions  and  heave  and  bow  accelerations  were  made, 
using  the  computer  program  for  comparison  with  the  experimental  results  of  Frldsma,^ 

Fridsma  tested  a series  of  constant-deadrlse  models  of  various  lengths  In  regular  waves  to 
define  the  effects  of  deadrlse,  trim,  loading,  speed,  length*to*beam  ratio  and  wave  proportions 
on  the  added  resistance,  heave  and  pitch  motions,  and  impact  accelerations  at  the  bow  and 
center  of  gravity.  Figure  3 shows  the  lines  of  the  prismatic  models.  The  models  were  towed 
at  CG  with  a system  that  permitted  freedom  in  surge.  The  computer  program  simulates  the 
model  being  towed  at  constant  speed  with  CG  at  the  baseline. 

Table  I presents  some  characteristics  of  the  model  and  experimental  conditions  for 
which  comparisons  were  made.  Most  of  the  comparisons  have  been  made  at  a speed-to-length 
ratio  V/y/'l  of  6.0  where  the  mathematical  model  is  expected  to  be  most  representative.  A 
limited  comparison  has  also  been  made  at  V/.y/T  ■ 4.0;  however,  no  comparison  has  been 
made  at  V/^C  ■ 2.0.  At  this  speed,  the  model  (or  craft)  operates  in  the  displacement  mode 
for  which  the  mathematical  formulation  is  not  valid. 

The  average  computer  run  corresponded  to  lO-sccond,  real-time,  model  scale;  however, 
only  the  last  2 seconds  were  considered  free  of  transient  effects.  An  example  of  the  compu- 
ter time  histories  of  pitch  and  heave  motions  is  shown  in  Figure  4.  Although  the  motions 
are  periodic,  they  are  not  perfectly  sinusoidal;  consequently,  in  determining  phase  relationship, 
the  peak,  positive-pitch  value  (bow  up)  and  the  peak,  negative-heave  value  (maximum  upward 
position  of  CG)  were  used  as  reference  points.  There  was  a difference  when  the  opposite, 
peaks  were  used. 
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TABLE  I - MODEL  CHARACTERISTICS  AND  WAVE 
CONDITIONS  FOR  COMPUTATIONS 

(Model  Length  ■ 1 14.3  cm  (3,75  ft);  L/b  ■ S:  - 0,608) 


CONFIGURATIONS 

SYMBOL 

wm 

LCG 

percent  L 

Radius  of 
Gyration 
percent  L 

Vly/T 

A 

20 

59.0 

25.1 

4,0 

B 

20 

62.0 

25.5 

6.0 

J 

10 

68.0 

26.2 

6.0 

M 

30 

60.5 

24.8 

6.0 

WAVE  CONDITIONS  FOR  CONFIGURATION  -- 

A 

B 

J 

M 

H/b 

IJL 

H/b 

\IL 

H/b 

III 

H/b 

0.1 11 

1.0 

0.1  II 

1.0 

0.111 

1.0 

0.111 

1.0 

0.111 

1.5 

0.111 

1.5 

0.111 

I.S 

O.III 

1.5 

0.111 

2.0 

0.1  II 

2.0 

0.1  II 

2.0 

O.III 

2.0 

0.1  II 

3.0 

O.til 

3.0 

O.III 

3.0 

O.III 

3.0 

0.1 11 

4,0 

0,111 

4.0 

0.1 11 

4.0 

O.III 

4.0 

0.1 11 

6.0 

0.222 

6.0 

0.111 

6.0 

0.111 

6.0 

0.334 

4.0 

0.111 

6.0 

J 

Corresponding  time  histories  of  bow  and  CC  accelerations  are  shown  in  Figure  5.  The 
bow  acccicration  was  computed  at  Station  0.  As  can  be  seen  in  these  piots,  the  impact 
accelerations  ranged  in  magnitude  from  cycle  to  cycle.  The  maximum  impact  for  negative 
value)  acceleration  computed  during  the  final  2 seconds  of  run  was  used  in  the  comparisons 
with  experimental  values.  In  some  instances,  particularly  near  resonance,  the  maximum 
impact  acceleration  was  more  than  twice  the  average  impact  value. 

Figure  6 shows  a comparison  of  variation  of  computed  and  experimental  pitch  and  heave 
motion  with  wave  height  for  the  20-degree  deadrise  model  In  a I S-foot  wavelength  and  fora 
speed-to-length  ratio  of  6.0.  Figure  7 shows  the  corresponding  impact  acceleration  at  the 
bow  and  CG.  The  computed  results  closely  follow  the  experimental  data,  except  for  CG 
acceleration  at  the  extreme  wave  height  condition,  where  the  computed  value  is  apparently 
much  lower.  Experimental  data  show  that  the  model  was  leaving  the  water  at  this  wave- 
height  condition.  The  computer  model  did  not  leave  the  water  but  came  very  close; 
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see  Figure  8.  Figure  8 is  a trajectory  of  the  computer  model  relative  to  the  wave  for  a 
selected  cycle  of  motion.  The  com(Aiter  model  behaves  very  much  as  expected.  On  the  left- 
hand  side  of  the  figure,  the  craft  is  planing  down  the  crest  ol  the  wave  and,  as  it  approaches 
the  wave  trough,  comes  very  close  to  leaving  the  water  before  slamming  and  submerging 
itself  deeply  into  the  front  of  the  oncoming  wave  crest. 

Figures  9 through  14  show  comparisons  of  the  computed  and  experimental  pitch  and 
heave  motions  at  \I\/T,  » 6.0  through  a range  of  wavelengths  and  at  a constant  wave  height 
of  2,54  centimeters  (I  inch)  for  deadrise  models  with  10.  20,  and  30  degrees.  The  data  have 
been  plotted  with  respect  to  the  coefficient  defined  by  Fridsma  as  L/X  (C^/(L/2b)^l 
Note  that  in  our  notation,  b is  the  half-beam. 

Comparisons  of  heave  and  pitch  for  the  10-degree  deadrise  model  shown  in  Figures  9 
and  10,  respectively,  show  excellent  results.  T he  computer  model  accurately  predicts  the 
secondary  peaks  in  the  pitch  and  heave  responses  at  0.19.  At  this  condition,  the  physical 
experimental  model  rebounds  so  as  to  fly  over  alternate  waves.  The  computer  model  oscillates 
at  half  the  wave-encounter  frequency  and  comes  close  to  leaving  the  water  at  alternate 
encounters  with  the  wave.  It  does  not  quite  leave  the  water  to  fly  over  alternate  wave  crests; 
nonetheless,  it  is  a good  representation  of  the  actual  motion. 

The  heave  and  pitch  comparison  for  the  20-degree  deadrise  model  at  V/^/T  " o.O  is  also 
excellent  us  cun  be  seen  in  Figures  1 1 and  12.  respectively.  No  experimental  phase  data  for 
the  condition  were  reported  for  Cy^  greater  than  0,072;  however,  extrapolated  results  (not 
shown)  are  in  line  with  the  computed  results,  The  pitch  and  heave  results  shown  in  Figures 
13  and  14  for  the  30-degrce  deadrise  model  are  good;  however,  responses  at  ■ 0.048  and 
» 0.072  are  higher  than  the  experimental  results. 

For  practical  considerations  a computational  scheme  for  planing  bout  motions  should  be 
valid  for  a range  from  approximately  V/>/T  ■ 4,0  to  V/>/T  ■ 6.0.  Computations  of  the 
motions  were  made  for  V/^/T  * 4.0  for  the  20-degree  deadrise  model;  see  Figures  15  and  16. 
Again  the  comparison  of  the  computed  heave  and  pitch  response  with  experimental  results  is 
excellent. 

Comparisons  of  the  computed  and  experimental  impact  accelerations  (or  largest  negative 
values)  are  presented  In  Figures  17  through  20.  Figures  17  and  18  show  bow  and  CG 
accelerations  for  the  10-dcgrcc  deadrise  model;  Figure  19  shows  similar  results  for  the  20- 
dcgrec  deadrise  model.  Figure  20  shows  the  results  for  the  30-degree  deadrise  model.  In  all 
cases,  the  comparison  appears  to  be  fair  to  good.  In  the  shorter  wavelengths,  X/L  = 1.0  and 
X/L  “1.5,  the  computed  accelerations  are  higher  than  the  corresponding  experimental  values. 
This  is  most  pronounced  for  the  1 0-degrec  deadrise  angle  model. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

A mathematical  model  of  a craft  having  a constant  deadrise  angle,  planing  in  regular 
waves,  has  been  formulated  using  a modified  low-aspect-ratio  or  strip  theory.  It  was  assumed 
that  the  wavelengths  were  long  in  comparison  to  the  craft  length  and  that  the  wave  slopes 
were  small.  The  coefficients  in  the  equations  of  motion  were  determined  by  a combination 
of  theoretical  and  empirical  relationships. 

A simplified  version  for  the  case  of  a craft  or  model  being  towed  at  constant  speed  was 
programed  for  computations  on  a digital  computer,  and  the  results  were  compared  with 
existing  experimental  data. 

The  comparison  of  the  computed  pitch  and  heave  motions  and  phase  angles  with  the 
corresponding  experimental  data  gave  remarkably  satisfying  results.  Comparison  of  the  bow 
and  CG  accelerations  was  fair  to  good. 

In  summary,  the  previously  described  mathematical  model  appears  to  be  a valid  represen- 
tation of  a planing  craft  in  waves  for  the  specific  craft  geometry  and  wave  conditions 
considered. 

To  make  the  computer  program  more  valuable  to  the  designer  the  following  additional 
work  is  recommended: 

1.  Improve  estimates  of  hydrodynamic  coefficients  to  obtain  better  acceleration  data 
and  to  include  more  complicated  ship  geometry. 

2.  Determine  added  resistance  In  waves, 

3.  Include  freedom  to  surge  and  to  add  components  of  propulsion. 

4.  Extend  to  the  case  of  Irregular  waves. 
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Figure  2 — Types  of  Two-DImensionnI  Flow 
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Figure  13  - Heiive  Response  for  30-Dcgrce  Deodrise  Model  nt  6.0 


Figure  14  - Pitch  Re^pon■e  for  30-Degree  Deadriie  Model  at  V/>/T>  6.0 
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APPENDIX  A 

EVALUATION  OF  HYDRODYNAMIC  FORCE  AND  MOMENT  INTEGRALS 

The  hyHrodvnamic  force  the  craft  experiences  in  the  vertical  direction  as  derived  in  the 
text  is: 


m,V-U  +rhgV  + CppbV2} 


cos  Q 


d{  + / apgAdf 
-'8 


where  U “ cos  0 - (/,-w^)  sin  d 
and 

V ® X(,j;  sin  0 + (z-w^)  cos  0 ~0( 

Another  force  acting  in  the  vertical  direction  is  the  weight  of  the  craft. 

The  first  two  terms  of  the  integral  are  evaluated  by  making  the  substitutions 


V ■ Xj;|j  sin  (/  - ii  J cos  0 - iVj,  cos  9 
+ 0(Xf.Q  cos  0 - sin  0)  + Wj,0  sin  9 


at' 

3? 


hi 

3« 


sin  0 


dw.  3w, 

— Lnw  - U — 

dt  '■  ^ 9{ 


and  noting  that 


3 m 

UV  “dt«-UVm 


auv 

3{ 


dj 


Using  the  previously  described  substitutions,  the  force  becomes 


■ j-Mg  cos  0 - Mg  sin  d icp^  +0,0  + MgO(z^,Q  sin  6 - cos  d) 

^ J vc^^  cos  0 (It  - ^fHg  WgS  tin  B d{ 

r ^'*'2  r 

■ / ~%r  0 dt  + / m,U  tt*  cos  8 dt 

yj  ■8?  •'{  ®« 


+ / apgAdt 
•^2 


This  is  essentially  the  form  in  which  the  integrals  have  been  computed  in  the  program. 

The  rate  of  change  of  the  sectional  added  mass  in  the  third  term  of  the  Integral 
expression  is  derived  by  relating  it  to  the  rate  of  change  of  depth  of  fluid  penetration  of  the 
section.  The  added  mass  of  a section  is  assumed  to  be  equal  to 


m,  - kg  it/2  pb^ 


for  which  the  time  derivative  is 


trig  ■ kgtrpbb 

where  b is  the  Instantaneous  half-beam  of  the  section,  and  kg  is  an  added-mass  coefficient, 
assumed  to  be  constant.  A value  of  k,  * i.O  was  used  in  the  computations  contained  in  this 
report.  For  sections  with  constant  deadrise.  which  is  an  imposed  limitation  of  this  work, 
the  half-beam  is  related  tu  the  depth  of  penetration  by 

b ■ d CO*  0 


.■^6 


where  d is  ilepth  of  penetration,  and  is  deadrise  angle. 

Taking  into  account  the  effect  of  water  pileup,  the  effective  depth  of  penetration  is, 
according  to  Wagner 


dg  = ff/2  d 


and 


b * dg  cot  (3  " ff/2  d cot  U 


where  7r/2  is  the  factor  by  which  the  wedge  immersion  Is  increased  by  the  pileup.  Using  this 
expression  for  the  half-beam,  the  rate  of  change  of  sectional  added  mass  becomes 


nij  “ kan'pb(tr/2  cot/3)d 

This  expression  is  valid  for  penetration  of  the  section  up  to  the  chine, 
exceeds  the  chine,  the  sectional  added  mass  is  assumed  to  be  constant, 


m,  - k 7r/2  pb^g^ 


When  the  immersion 
i.e.. 


where  is  the  half-beam  at  chine. 

The  submergence  of  a section  in  terms  of  the  motions  is  given  by 

h = / - r 


where  /.  = { sin  ^ + f cos  U 

r = r^i  cos  ^ cos  -i-  f '-in  d)  + cot } 

For  wavelengths  which  are  long  in  comparison  to  the  draft  and  for  small  wave  slopes,  the 
immersion  of  a section  measnreil  perpendicular  to  the  baseline  is  approximately 

/ - r 

d 

cos  0 - u sin  0 


where  v - wave  slope 


The  rate  change  of  submergence  d is  given  by 


^ ^ i - r ^ (z  - r)  _ d(cosg  - using) 

"cos0-usin6  (cose -u  sin  0)2  ' 

Since  immersion  (z-r)  Is  always  small  in  the  valid  range  of  the  previously  described 
expression,  the  relationship  can  be  further  simplified  to 

— LJ. — 

cos  e “ u sin  6 


and 


nig  ^ k,ff/)b(jr/2  cot/J) 


(i  - h 

cos  e - u sin  e 


The  expansion  of  the  integral  expression  for  the  hydrodynamic  moment  In  pitch 
follows  the  procedure  used  for  the  vertical  force.  The  results  are  summarized  as  follows 

Fg  • - i,e  + Qg  cos  e - Qgflfzcc  e - x^c  cos  e) 
dw 

sin  e w^fdf 


+ f Vrh,«d{+  /pCpbV2{dt 
-'fi  -'fi 

+ m,UV{|  + / mgVUdt 

® 'stern  ® 

r 3w, 

Vj  IT  “ 

r 


cos  e {dt 


+ / apgA  cos  0 td{ 

A 


The  only  additional  moments  are  the  buoyancy  moments.  All  other  moments  are  considered 
to  be  zero  for  the  specific  problem  considered  in  this  report. 
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APPENDIX  B 

COMPUTER  PROGRAM  DESCRIPTIONS 


OVERVIEW 

The  equations  of  motions  developed  In  the  previous  sections  of  this  report  have  been 
solved  by  means  of  digital  computer  programs.  Two  major  programs  have  been  developed: 
the  first  (MAIN)  solves  the  equations  of  motion  using  the  Runge-Kutta-Merson  integration 
algorithm  and  generates  time  histories  that  are  stored  on  the  system  disk.  The  second 
(PLTHSP)  generates  California  Computer  Products  Company  (CALCOMP)  pen  plots  from  the 
disk  files.  All  programs  were  designed  to  operate  on  the  Control  Data  Corporation  computer 
system,  located  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  in 
Carderock,  Md. 

Descriptions  of  input  data  required  to  execute  the  programs,  job  control  cards,  and 
programs  follow.  Sufficient  detail  is  presented  for  this  appendix  to  serve  as  a manual  for 
use  and  maintenance. 


JOB  CONTROL  CARDS  FOR  PROGRAM  MAIN 

Job  control  cards  for  program  MAIN  which  computes  time  histories  of  the  motion 
variables,  are  described  as  follows.  If  CALCOMP  plots  are  not  desired,  TAPES  need  not  be 
cataloged. 

Job  Control  Language  Card: 

Job  Card 
Charge  Card 


REQUEST, TAPE9,*PF. 

REQUEST,TAPE2,*PF. 

REQUEST  ,TAPE4,*PF. 

ATTACH.BINAR,SEFZARNICKNEWB. 

ID-XXXX. 

ATTACH,NSRDC. 

LDSET(LIB-NSRDC). 

BINAR. 

REWIND,TAPE2. 

REWIND,TAPE4. 

COP  Y(TAPE2, OUTPUT) 

COPY(TAPE4, OUTPUT) 


Comment 

Standard  facility  card 

Standard  facility  card 

Reserves  space  for  CALCOMP  plot  dutu 

Print  output  file  1 request 

Print  output  flic  2 request 

Attaches  binary  run  file 

Attaches  library  routines 
Loads  library  routines 
Loads  und  executes  run  file 

Rewinds  time-history  files  for  printing 

Prints  tinie-hlstory  file 
Prints  time-history  file 


Comment 


Job  Control  Language  Card: 

CATALOG.TAPE9,  SEFZARNICKDATA. . , Catalogues  file  for  plot. 

ID-XXXX.  (SEFZARNICKDATA  CAN  BE  ANY  NAME) 

7/8/9  END  OF  RECORD 

DATA  CARDS  (1-5) 

6/7/8/9  END  OF  FILE 

INPUT  DATA  CARDS  FOR  PROCPAM  MAIN 

Input  data  used  by  program  MAIN  are  read  from  data  cards  in  NAMELIST  and  In 
standard  format.  A descriptor  c'  the  FORTRAN  symbols  appearing  in  NAMELIST  follows. 
For  simplicity  in  the  text  that  follows,  It  is  assumed  that  NAMELIST  input  occupies  only 
one  card.  More  cards  can  be  used  if  necessary. 

Card  KNAMELIST  FORMAT,  / / 1 

A The  absolute  error  for  KUTMER  (six  values) 

NPRINT  lf"l,  print  normal  output 

If" 2,  matrix,  inverse  matrix,  F-column  matrix,  and  KUTMER  results 
lf"3,  integral  results 

If "4,  calculated  values  constant  for  given  input  values 
NPLOT  lf"0,  no  plot 

If" I,  printer  plot  of  results 
END  Number  of  runs  to  be  made 

W Weight  of  craft  In  pounds 

BL  Boat  length  in  feet 

TZ  Thrust  component  In  z direction 

TX  Thrust  component  in  x direction 

XECG  Distance  from  center  of  gravity  to  center  of  pressure  for  drag  force  in  feet 

XP  Moment  arm  of  propeller  thrust 

XD  Distance  from  center  of  gravity  to  center 

DRAG  Friction  for  drag  force 

RO  Wave  height 

LAMBDA  Wavelength 

RG  Radius  of  gyration  In  feet 

T Propeller  thrust  In  pounds 

GAMMA  Propeller  thrust  angle  in  degiees 


Card  1 (Lontinucd) 


ECG 

Longitudinal  center  of  gravity 

NCG 

Vertical  center  of  gravity,  nondimensionalized  by  ship  length 

KAR 

Added-mass  coefficient 

BETA(I) 

Dead-rise  angle  in  degrees 

ESTd) 

Station  position  in  feet 

NUM 

Number  of  stations 

XA 

Initial  time 

XE 

Stop  time 

HMIN 

Minimum  step  size 

HMAX 

Maximum  step  size 

EPS 

Error  criterion 

Card  2 (Format  8F10.0) 

(X(I), 1-1.6) 

Initial  conditions 

X(l) 

Velocity 

X(2) 

Z 

X(3) 

0 

X(4) 

X 

X(5) 

z 

X(6) 

0 degrees 

Card  3 (8F10.0) 

START 

Time  to  turn  on  (RMP)  function  (see  page  48) 

Risi; 

Duration  of  RMP 

Card  4 (BF10.0) 

TME 

Time  at  which  integration  interval  is  to  he  chonged* 

MMX 

New  maximum  interval  size  after  TMF. 

HMN 

New  minimum  interval  size  for  KUTMER  to  subdivide 

♦If  tlili  option  l5  not  uied  lel  TMli  to  »top  timo  on  run. 
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Card  5 (8F10.0) 


PERCNT  Percentage  of  bout  length  subtracted  from  longitudinal  center  of  gravity  to 
obtain  X - point  where  acceleration  computations  are  made 


JOB  CONTROL  CARDS  FOR  PROGRAM  PLTHSP 

Job  control  cards  for  program  PLTHSP  which  generates  CALCOMP  plots  of  time 
histories  computed  by  program  MAIN  are  described  in  this  section. 


Job  Control  Language  Card: 

Job  Card 
Charge  Card 
REQUEST, TAPE7.HI. 
VSN(TAPE7-CK0323). 

ATTACH.CALC936. 

attach.binar.sffzarnickplotb, 

ID-XXXX. 

LDSET(LIB-CALC936) 

BINAR, 

7/8/9  END  OF  RECORD 
DATA  CARDS 
6/7/8/9  END  OF  FILE 


CamjngnJ 

Standard  facility  card 

Standard  facility  card 

Tape  for  CALCOMP  plot  data 

Volume  serial  number  of  tape  for 
CALCOMP  plot 

Attaches  CALCOMP  library  routine 
Attaches  plot  program  run  file 

Loads  CALCOMP  library  routines 
Runs  plot  program 


INPUT  DATA  CARDS  FOR  PROGRAM  PLTHSP 

Two  or  three  data  cards  are  made  ready  by  PLTHSP,  depending  on  the  options  selected. 
Standard  input  format  is  employed.  A description  of  the  necessary  data  cards  follows. 


Card  1 (8F10.0  Format) 

XAXIS  Lengtli  of  x axis  in  inches 

YAXISP  Height  of  pitch  component  axis  in  inches 

YAXISH  Height  of  heave  component  axis  In  Inches 

HT  Height  of  lettering  in  inches 


Card  2 (110 

lA 


Format) 

lf*0,  no  plots  for  bow  acceleration  and  center  of  gravity  acceleration 
If“l,  plots  previously  mentioned  information 
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Card  3 (8F10.0  Format)  - Only  Necauary  If  lA*  1. 

YAXISB  Height  of  bow  acceleration  axis  in  inches 
YAXISC  Height  of  CG  acceleration  axis  in  inches 

PROGRAM  MAIN 

Program  MAIN  reads  all  necessary  input  data  from  cards,  sets  up  initial  values,  computes 
constants,  calls  KUTMER  to  determine  the  state  variables  at  TIME  for  the  period  from  XA 
to  XE  In  increments  of  HMAX.  A table  state  variables  is  created  for  every  PTIME-th  value. 
The  values  for  X/H  and  0p/2ffHA  are  calculated  and  printed.  If  the  plot  option  is  on.  a 
printer  plot  will  be  produced. 

Subroutin*  COMPUTIX) 

This  routine  computes  pitch  moment  NL  and  lift  force  FL.  excluding  added  muss  terms, 
using  values  of  integrals  computed  in  subroutine  FUNCT,  The  argument  X contains  the  state 
vector. 

Subroutine  DAUX 

This  subroutine  Is  culled  from  KUTMER  or  EULER.  It  determines  the  values  of  m^.  b. 
and  bl*,  based  on  the  following  equations 

h,^fn  ■ -{(I)  sin  0 + ffh  cosO  - rd) 

where  r(l)  ■ r„  cos  k + {(1)  cos  0 + f(l)  sin  0 + ct) 

Then  for 


h^tlt  > 0. 

lUtll 

dtl)  • — 

cos  0 - il)  sin  0 

where  V(l)  “ - r„k  sin  0 (x^.^.  + ^(  1)  cos  d + (1)  sin  0 + ct) 

If 

dtll  > tun  (dd)  2/7r} 
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set 


nijd)  - 

b(l)  - b^(l) 

bid)  “0 

- k(l)(p/2)nb*  (I) 


d(I)<b_(I)  tan  W(I))(2/jr) 


bd)  -d{t)cot(^(i))(7r/2) 
bid)  - b(l) 

m,d)  ■ k,(l)(p/2)irb^(l) 


h^dXO; 

nigd)  ■ 0,  b(I)  ■ 0,  bid)  ■ 0 


Thli  lubroutlne  then  calli  FUNCT  which  in  turn  calls  COMPUT  to  determine  the  values 
of  and  the  lift  force  and  moment.  The  values  of  and  F^  are  used  to  compute 
the  following 


F,  ■ + F^  sin  0 - D cos  9 

Fj  ■ Tj  + Fj^  cos  0 + D sin  0 + W 

• NL-Dxd+T,p 


*bt  arrsy  Is  set  up  for  Intagratlnni  for  portion  of  hull  for  which  chine  li  not  Immaried, 
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The  mass  inertia  matrix  is 


A 

A 

A 

A 

A 

A 

A 

A 

A 


11 

12 

13 

21 

22 

23 

31 

32 

33 


M + M,  sin^  S 
Mg  sin  6 cos  0 
-Qg  sin  d 

'^12 

M + Mg  cos*  6 
- Qg  cos  8 


The  matrix  is  inverted  by  the  system  routine  MATINS,  The  Inverted  matrix  is  then 
used  to  solve  the  following  equations  which  determine  the  state  vectors. 


^co  * ^11*  ^1  ^ '^la**  ^2  ^13*'  ^3 

8 “ Aj,*‘ F,  + Ajj  * Fj  + Aj3*‘ Fj 


Subroutine  FUNCT  (X) 

This  routine  evaluates  various  integrals  appearing  in  the  force  and  moment  inathcmuticul 
models.  The  integrals  are  evaluated,  using  a trapezoidal  integration  algorithm.  The  argument 
X contains  the  state  vector.  A list  of  Integrals  that  are  evaluated  is  presented. 
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Subroutin*  INPUT 

This  routine  reads  In  NAMELIST/HSP/  which  contains  the  Initial  data  concerning  the 
craft  and  sea  conditions  pertinent  to  all  the  runs  to  be  made,  it  is  set  up  so  that  most  of 
the  data  are  given  default  values  by  means  of  data  statements  in  subroutine  INPUT.  These 
data  statements  can  be  overridden  during  execution  by  reading  values  In  on  cards.  For  further 
explanation  of  the  specific  variables  see  section  on  the  input  data  cards. 

This  routine  also  "initializes"  constant  such  as  w.  p,  and  g.  It  uses  the  input  values  to 
calculate  the  keel  profile  and  planform  arrays,  NO  and  BM,  wave  constants,  system  mass  and 
inertia,  and  maximum  mass  and  depth  of  chine  at  each  station. 

Subroutine  KUTMER  (NEQS,  TIME,  HMAX,  X,  EPSE,  A,  HMIN,  FIRST) 

This  is  a Runge-Kutta-Merson  Integration  routine  that  is  capable  of  changing  the  si/e  of 
the  interval  over  which  it  integrates  to  meet  specified  error  criteria.  It  is  therefore  an 
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accuia’c  method  for  a system  that  may  oscillate  more  rapidly  than  the  initial  integration 
intcrviil,  A minimum  step  size  prevents  the  routine  from  subdividing  the  interval  indefinitely. 
The  input  arguments  are: 

NEQS  Number  of  dependent  variables  in  the  x array 

TIME  Actual  time  (independent  variable) 

HMAX  Increment  for  which  the  solution  is  to  b”  returned 

X Vector  of  dependent  variables 

EP3E  Relative  error  criteria  specified  for  each  component  of  x and  used  for  the 

components  of  x less  than  the  absolute  value  of  A 

A Absolute  error  criteria 

HMIN  Minimum  step  size  allowed 

FIRST  Set  to  zero  on  first  call;  a value  of  1 is  assigned  by  KUTMER  on  subscijuent 

calls  tor  which  the  error  criteria  are  satisfied,  otherwise  a value  of  2 is 
assigned 

Subruutin*  PLOT2  (F,  FMIN,  FMAX,  NVAR,  NFUN,  N1,  N,  XC,  DELX) 

Data  stored  in  the  two-dimensional  array  F are  plotted,  using  the  printer  by  subroutine 

PL0T2.  As  many  as  2C  different  functions,  having  evenly  spaced  abscissa  values,  can  be 

plotted,  The  output  is  written  on  Unit  6.  A des'iript'on  of  variables  follows. 

F Array  containing  data  to  be  plotted;  the  Jth  point  of  the  1th  function  is 

stored  in  F(I,J) 

FMIN  An  array  of  minimum  lunctional  values;  the  minimum  of  tlie  Ith  function 

li.  stored  in  FMIN(l) 

FMAX  Same  as  FMIN  only  for  maximum  values 

NVAR  An  array  of  titles  for  each  function  to  be  plotted 

NFUN  Number  of  functions  to  be  plotted 

N1  First  dimension  of  array  F 

N Number  of  points  to  be  plotted 

XO  First  abscissa  value 

DELX  Abscissa  increment 

Subroutine  PLOTER  (FX,  XA,  HMAX,  LAMBDA,  IB,  NWAVE) 

The  routine  initializes  various  values  required  to  generate  printer  plots  and  computes 
pitch-aiid-heave  ratios.  The  printer  plots  that  arc  generated  consists  of  pitch-and-heave  time 
histories.  A description  of  input  variables  follows. 


FX 

XA 

HMAX 

LAMBDA 

IB 

NWAVE 


A two-dimensional  array,  containing  time  histories  to  be  plotted 
Initial  time 

Time-interval  ineioment;  time  interval  between  values  in  FX  is  given  by 
HMAX*PTIME 

Wavelength 

Nu, Tiber  of  values  to  be  plotted 

Position  in  FX  at  which  wave  is  completely  turned  on 


Function  RMP  (T.  START.  RISE) 

The  RMP  is  a function  that  calculates  a value  oetween  0 and  I corresponding  to  time  T, 
based  on  a straight  line  from  time  START  with  a vaiue  of  0 to  time  START  plus  RISE  with 
a value  of  i.  It  Is  used  to  lower  the  initial  wave  amplitude  to  avoid  large  transients  at  start 
of  llii  <'omputai:ions. 

The  arguments  are; 

T Actual  time 

START  Time  at  which  to  begin  the  ramp  from  0 to  1 
RISE  Duration  of  rise  from  0 to  I 

The  function  reaches  the  value  1 at  time  START  plus  RISE,  If  the  rise  Is  0.0,  RMP  will 
return  a value  of  O.S. 


Subroutine  TRAP  (F,  DX.  NPTS,  ANS) 

This  routine  performs  the  evaluation  of  an  Integral  using  a trapezoidal  approximation. 
The  argument  variables  are  defined  as  follows: 

F Array  of  Int  egrand  values 

DX  Increments  at  which  F is  evaluated 

NPTS  Number  of  values  in  F 

ANS  Result,  which  is  equal  to 


DX 


F(i)-0,5(F(I)  + F(NPTS)] 


PROGRAM  PLTHSP 

This  program  uses  a data  file  created  by  program  MAIN  to  create  CALCOMP  plots.  The 
data  are  read  from  logical  Unit  9 and  are  rewritten  on  Unit  7 for  CALCOMP  input.  Program 
PLTHSP  sets  the  tape  output  unit  equal  to  7 and  calls  SUBROUTINE  CALPHI  to  execute  the 
plot  procedures. 
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Subroutine  CALPLT 

This  subroutine  manages  all  the  I/O  operations  and  performs  the  necessaiy  calculations 
required  to  generate  the  plots.  After  reading  the  card  data  (two  or  three  cards)  subroutine 
READT  is  called  to  read  the  data  file  (Tape  9)  created  by  program  MAIN.  The  CALCOMP 
initializing  routines  are  called  next,  after  which  a call  to  subroutine  ESCALE  calculates  the 
necessary  scaling  factors.  Subroutine  EXAXIS  is  called  next  to  determine  the  placement  of 
the  plot  tick  marks  and  identifying  digits.  The  CALCOMP  plot*generation  subroutines  are 
now  called  and,  depending  on  the  option  ocflned  by  the  lA  parameter  on  card  2,  plots  of 
pitch  and  heave  at  the  bow  and  CG  location  are  generated  as  functions  of  time  if  lA  ■ 1. 

Subroutine  EAXI8 

The  subroutine  is  analogous  to  the  CALCOMP  AXIS  routine,  The  only  exception  is  that 
the  tick  marks  are  not  necessarily  inch,  and  the  height  of  the  characters  is  defined  by  the 
input  parameter  HT,  Function  NDIGIT  is  called  to  determine  the  number  of  digits  necessary 
to  print  an  even  increment  of  the  plots  functions  on  the  axis. 

Subroutina  ESCALE,  ADJUST,  and  FUNCTION  UNIT 

These  subroutines  find  the  scale  to  be  used  on  the  plot  axis.  Function  UNIT  is  called 
to  determine  the  axis  Increment  size  after  which  subroutine  ADJUST  is  called  to  extend  the 
minimum  (AMIN)  and  maximum  (AMAX)  values  so  that  they  are  even  multiples  of  the  axis 
increments. 

FUNCTION  NDIQIT 

This  function  finds  the  number  of  digits  necessary  to  print  even  increments  of  the 
function  on  the  axis.  Both  the  number  of  nlaces  in  tho  entire  number  (NDIGIT)  and  tlie 
number  of  decimal  places  (ND)  are  determined,  after  which  l>ie  value  of  each  increment  on 
the  axis  (ANUM)  is  calculated. 

Subroutine  READT 

This  subroutine  reads  the  data  file  cieated  by  program  MAIN.  Data  file  records  are 
read  until  the  message  end  of  file  is  encountered.  Each  recoid  is  read  in  the  same  format  as 
it  was  written  in  MAIN.  The  information  is  printed  to  allow  the  user  to  inspect  the  created 
nie. 
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n r>  o oo 


BEST  AVAIWBIE  COPY 

LISTING  OF  COMPUTER  PROGRAM  FOR  MOTION  COMPUTATIONS 


PROORAM  MA IN ( INPUT fOUTPUT » T APt5«INPUT * TaPC6»0UTPUT t ’ APe3>Sl 2 1 

• TAPe2aS12<TAPEAaSI2*TAPE9> 

REAL  lT«K,LAHeOA,M,MA,MNAA«N.NCO.NU«MASS«NL«lA«KAR 
INTCaCR  ENO 

DIMENSION  X(6)trX(2>A00) 

COMMON  /CONST/  NCO«ECO«Pl«OPR«RPO»ORAVrY«RHO«K«NUMtMA(120) «C0»TA, 

• B(120).BETA|Hwa20)*T2«UHAG«N,xO«T«XP*M.lT« 

• DELTASiTXiESTIliOlfCtROtKARtMMAXIl  0)«TeSTU20)i 

« NIlEOliPMAL^ 

COMMON  /SHIP/  MASS«CINTiOA»CG«CC2*Ce3«OMUiCOMU«E20NU»C3DMUtBPiBMM 

• NLiPLiIAfElUOl 
COMMON  /IN/  BM(120) •811120) «VLLIN 
CUMMON/UUT /NPRl NT •NPLOT teNO 
C0MN0N/TrRMS/Tl«T2tr3«T4«TS,TA«T7.rB 
COMMON  /SCAHAVE/  START tRl SCtRAMP 
COMMON  /INTCH/  It iKTTIlO) tOIPP <10) 

COMMON  /IN2/  N0020)iXA«XC*HMAX«HM|N,A<6)»EPSe(S) •lambda 
COMMON  /ACCEL  / XACCL^BMACL^COACL^BL 

CALL  INPUT 

COMPUTE  INTEORATION  INTERVAL  INFORMATION 

NLCSS  ■ NUM«1 

I • I 

II  ■ I 

differ  ■ csT<i*ii-esT(n 
KTTIIll  • I 
OIFF(Il)  ■ differ 
DO  2S  I«2^NLESS 
OlFFERa  eST<l*l)-EST(I) 

XTTdl)  ■ KTT<ll)*l 
IFIOTrFtN.NE.UlFFdlDOO  TO 

00  ro  23 

2A  It  ■ IIM 
NTTdl)  ■ \ 

DlFFdl)  ■ differ 
a CONTINUE. 

KTTdl)  • KTTdI)*l 

• • • • • CHECK  IF  NUMBER  OF  INTERVALS  EXCEEDS  DIMENSION 

IF  dItOT.lOl  PRITE'6»2S)  (KTT ( 1 ) .OIFF ( 1) . I»l . 1 1 ) 

IFdl.OT.lO)  STOP  A 

• • • • • POINT  AT  XMICH  MULTIPLE  RUNS  START 
a continue 

TIME-XA' 

KOUNT*! 

ENO-END'l 
MRITE<6i39) 

39  format  dHI) 

••••••••  READ 

Xd)  • VELOCITV* 

X<A)  ■ X, 

THETA  IS  pEAU  in 


DOT 


IN  initial  conditions 

X(2)  a 2 OoT*  X<3)  ■ THETA 
X<8)  a 2*  XI6)  a theta 

oeorees  then  converted  TO  radians  in  program 


REAO<S»10) (X<I) •lal«6« 


main 
MAIN 
main 
main 
main 
main 
main 
main 

MAIN  10 
MAIN  II 
MAIN  12 
main  13 
•MAIN  lA 
main  15 
main  16 
main  17 
MAIN  IB 
MAIN  19 

main  20 
MAIN  21 
main  22 
MAIN  23 
main  2A 

MAIN  2S 
MAIN  26 
main  27 
MAIN  26 
MAIN  29 
MAIN  30 
MAIN  31 
MAIN  32 
MAIN  33 
MAIN  3A 
MAIN  3S 
MAIN  36 
MAIN  37 
MAIN  36 
MAIN  39 
MAIN  AO 
main  Al 
MAIN  A2 
MAIN  A3 
MAIN  AA 
MAIN  AS 
MAIN  A6 
MAIN  A7 
MAIN  A8 
MAIN  A9 
MAIN  50 
MAIN  51 
MAIN  52 
MAIN  53 
MAIN  5A 
MAIN  55 
MAIN  56 
MAIN  57 
MAIN  56 
MAIN  59 
MAIN  60 
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BEST  AVAIUBLE  COPY 


DATA  » useo  IN  RANP  FUNCTION*  TO  TURN  ON  WAWt  N*1N  61 

R£A0(5,I0»STAHT.«ISC 

MAiN  o3 

10  FORNATjBFH.A) 

• •••••••  HR  lie  OUT  TM€  INPUT  VAtUCS  **AtN  65 

HRlTt<6.l9)  START, RISC .KAR  ^ 

19  FORMAT  (•'  start  • WISE  ■ ",F10iA,/,“  RAH  • ''.FIONAIN  6T 

..4»  »• 

MAIN  69 

TME  Is  THE  TIME  AT  WHICH  THE  INTEGRATION  INTERVAL  ts  MAIN  TO 

TO  NC  changed 

HMA  Is  THE  NEW  MAKlMuH  INTERVAL  SUE  AFTER  TIME  THE  MAIN  TE 

HMN  Is  THE  NEW  MlNlNgN  INTERVAL  SUE  FOR  KUTMER  TO  SUB-OIVIOE  MAIN  T3 

THE  MAXIMUM  INTERVAL  uR  TO  MAIN  T4 

IF  TMIS  option  is  not  used  set  TME  TO  THE  STOP  TIME  OF  TME  RUN  MAIN  75 

READIS.IO)  TME, MMX, HMN  JT 

WRITE(6»IU  T<tE,HMAt,MMX,MM|N,MHN  MAIN  7# 

U FORMAT!*  AT  TIME  *,77,1, • THE  MAXIMUM  INTERVAL  SUE  FOR  INTEORATIMAIN  T9 

•ON  MILL  RE  CMANGCO  FROM  •,FI0*4,*  TO  •tFIO.A,/,  MAIN  SO 

••  AND  TME  MINIMUM  SUE  FOR  HALVING  CHANGES  FROM  •,F10,4,»  MAIN  Ol 

• • TO  *,FlO,4l 

ADJUST  TME  TIME  FOP  CHANGE  OF  INTEGRAUON  INTERVAL  MAIN  B3 

FOR  CHECK  AGAINST  TIME  IN  THE  INTEOARTION  LOOP  MAIN  S4 

TM  ■ TMC>(HMAX/t*>  MAIN  05 

SET  SWITCH  FOR  CALCULATION  OF  PITCH  AND  HEAVE  RATIOCS  MAIN  S6 

ON  NEXT  CALL  TO  PLOTER 

IPT  » 0 •• 

irCTME.EO.xC)  IPT  • I JJUJJ  •« 

READ(S,10>  PERCNT  Si 

XACCL  ■ ECf»-PERCNT*»L  ** 

WRITE(6«U)  PCrCNT, XACCL  93 

12  FORMAT!*  THE'X  USED  FOR  THE  BOv  ANO  CO  ACCELERATION  COMPUTATIONS  MAIN  94 

419  COUAL  TO  EC0-«*F19.A*TH*BL  OR  »FI0,4>  MAIN  95 


WR1TC!6«23) 

HR1TE(6«47) 

23  FORMAT !lH  ,//} 

47  FORMAT!"  STATION  NO," ,3X,"0EAO  Rl5C",8X,"C5T",8X,‘*N0", 
* I0X,"BEAm*'» 

WRITE (6»55)  1 UtBCTx.tSTdl  ,NO!I),BM!l))  ,I*I,NUM) 

55  FORMAT (6x,12,5x,Fl0.4,4X»FlB.4,4X»f I6,4,3X,ri0.4» 
WRITE(6f23) 

WRITE  (6,S6I  <X<I),1-I,6> 

56  FORMAT!"  X VALUES", 4X,6!FI0,4»2X»» 

••••••«•  CHXNOC  INPUT  From  OEGRELS  TO  RADIANS 

X(31  ■ X(3»*RP0 
X!6)  • X!6)*RP0 

have  • STAbT*RI$E 
NWAVe  a 0 

••••••••  write  OUT  COMPUTCO  ARHATS 

WRlTE(6,57iM,lT,K,C,PHALF ,P:,GRAVTT 
IF(NPRINT,LT.4)  go  to  62 
WRITE  !6,5A)  !e(I),I*I*NUH) 

WRITE  !6,S9)  (N< I t « I*l ,NUM) 

WRITE  (6,64)  (HMAX<1) flaltNUM) 

WRITE  (6»64)  ITCST(I).Ial«NUM> 


MAIN  95 
MAIN  96 
MAIN  9T 
MAIN  9B 
MAIN  99 

main  100 
main  101 
main  102 

MAIN  103 
MAIN  104 
MAIN  105 
MAIN  106 
MAIN  107 
MAIN  lOB 
MAIN  109 
MAIN  IIO 
MAIN  111 
MAIN  U2 
main  U3 
main  114 
MAIN  115 
MAIN  116 
main  117 
main  IIB 

MAIN  119 


fc*'-.;  ;• 


BBT'AVAIUBLE  copy 


continue 

MRITE(6«28)  <KTT<I>«01Pr(I)tIalin> 

IS  fORNAT(«  KTTtOIPf  tt IlOf 2X«r 10.4) 

Sr  ronMAT(4N  mb  «flO*A«4H  !■  *F10.«*4H  kb  «riO.A,4H  Ca  »ri0.4«UH 
«RHO/|a  tnOtAtSM  Pla  .PIO.AiIAh  ORAVlTYa  .riO.41 
$S  rONNAT  <•'  ta)“*10Pl0.A) 

$9  PORNAT  N<1)*«*10P10.A) 

6*  PONMAT  HNAX(t)'Sior|0,4) 

60  PORNAT  TC9T(l)HtlOP10.4) 

le  • 1 

IRRINT  • NPRINT 
WRITE(A«91) 

Cbbbbbbbb  mRITE  HEAOINOS  AND  CONDITIONS  AT  TINE  > 0. 

91  PORNAT (lHl,|X<‘>TlME»«9X*«X00T'S9X«»l00Ti*t9X»HTHCTA  o6T»>6X* 

« lHX«9X«lHZt9X»SHTHCTA«9Xi2HNL,9A«2HPL* 

« AXtSHUUW  ACCL.AXtTHCO  ACCLt//) 

WRITE  (A«92l  TlMEi  (X(I|  tlal  »6l  tNL*Pl.,»WACi.»CaACt 
WR1TE(9>  TINei(X(ll »laA»6) iSWACCtCOACC 
KOUNT  a K01)NT*1 
PXa«lB)aX(S) 

PX<2«IB>aX(6> 

XKUTMa(XC-XA>/HMAX*,OS 

IKUTN  B (TNC-XA)/HMAX  * (XE«THe>/NHX  « «OS 
PIRSTbB.O 
NEOSbO 
IKUTSaO 
C 

C start  op  INTEORATION  LOUP 

c 

BSl  continue 

NPRINT  • IPRINT 

CB«BBB«Bt  CHECK  pitch  .OT,  tSESB  RAOIANS 
IP(X<6).OT. *5136)00  TO  BS3 
CBBBBBtBB  PERPOPM  INTEORATIONS 

IPaiME.LT.TM.OR.TNE.EO.xC)  OU  TU  96 
IMIPT.EU.I)  GO  TO  9B 

HMIN  a hMN 
HNAX  a HNX 
FIRST  a B.O 
9B  CONTINUE 

CALL  KUTMEo<NEOS«TINCtHMAX*XtEPsE»A*HNlN*PlRST) 

IKUTSalHllTS*! 

IPIPIRST.E'J.DGU  TO  B61 
IP(K0UNT.NF.l.AN0.K0UNT.Ne.4|)  CU  TO  99 
HRITE<4t9)) 

KOUNTb) 

Caaaaaaa*  WRITE  OUT  TINE  INTERVAL  RESULTS 
99  WRITE(4i92)  TINE. (X (I ) , lal ,6) tNLtrL*BWACL*CCiACL 
WRITE <6.93) TliTE*T3»T4*T9tT6.TT,TB.BMM,WP 
WRITE  (9)  TINE.CXd)  .Ia4*6)  tBNACL.COACL 
IP<TIMC.LT.TH.OR.TMe.EO*xe)  OU  TO  100 
IP(IPT.EO.I)  Ou  TO  200 

CALL  PLUTC’lPX.XA.HNAX.LANeOA.lB.NWAVE.lPT) 

IPT  a 1 

19  a 0 

XA  a TIME 

first  a 0.0 

HMIN  ■ HMN 
HMAX  a MMX 


HA 

IN 

120 

HA 

IN 

121 

HA 

IN 

122 

ha 

IN 

123 

HA 

IN 

124 

ha 

IN 

125 

HA 

IN 

126 

HA 

IN 

127 

HA 

IN 

12B 

ha 

IN 

129 

HA 

IN 

130 

HA 

IN 

131 

HA 

IN 

132 

HA 

IN 

133 

HA 

IN 

134 

HA 

IN 

135 

HA 

IN 

136 

HA 

IN 

137 

HA 

IN 

13B 

MA 

IN 

139 

HA 

IN 

140 

HA 

IN 

141 

MA 

IN 

142 

MA 

IN 

143 

HA 

IN 

144 

ha 

IN 

145 

ha 

IN 

146 

HA 

IN 

147 

HA 

IN 

14B 

HA 

IN 

149 

ma 

IN 

150 

HA 

IN 

151 

HA 

IN 

152 

MA 

IN 

153 

MA 

IN 

154 

MA 

IN 

ISS 

MA 

IN 

156 

HA 

IN 

1S7 

MA 

IN 

19B 

HA 

IN 

159 

ma 

IN 

160 

MA 

IN 

161 

MA 

IN 

162 

MA 

IN 

163 

MA 

IN 

164 

MA 

IN 

16$ 

MA 

IN 

166 

MA 

IN 

167 

MA 

IN 

16B 

MA 

IN 

169 

HA 

IN 

170 

MA 

IN 

171 

MA 

IN 

172 

MA 

IN 

173 

MA 

IN 

174 

MA 

IN 

175 

MA 

IN 

176 

MA 

IN 

177 

MA 

IN 

178 
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nooririonrtonoorto 


BEST  AVAIUBIE 


{fO  CONTINUE 

iB"ia*i 

rx(CtiB)*x(6) 

93  POBNXTC  •MOeiO.Ai) 

92  rOBMATIlX»U(^10.<>«2X)) 

I BO  CONTINUE 
‘ KOUNTbKOUNT*! 
ir(NNAVE,OT.B)00  TO  21 
ir<TlHE.OT,WAve>NMAVC«KOUNT 
21  continue 

irETIHE.LE.xe. AND, IKUTS.lt. 1KUTM)0U  TO  BSl 
MRlTe<2»aS?> 

BBA  CONTINUE 

BBB  rORMAT(‘»  CNO  OF  KUTNCR**) 

BB3  CONTINUE 

CALL  PLUTE*)  (FX « XA .HMAX .LANBOA » IB • NM AVE • IPT  > 

CHECK  for  LAST  RUN  IF  NUT  CYCLE  HACK  TO  READ 
C NCR  DATA  FOR  NEXT  RUN 
IFEENO.NE.DOU  TO  B 
SO  TO  999 

C • • • • KUTHER  ERROR  NESSAOES 
S«l  WRITE<0«S«>) 

8B2  FORMAT ("  ERROR  CRITERION  IN  KUTMER  CAN  NOT  BE  MET*') 

WRITE  (BiSA)  (X(nil>t*BI 
WRITE  («iSA»  TIME 
•6  FORMAT  <*•  TIME  ■•••FIO.A) 

IFEENO.NE.UOb  TO  B 
00  TO  003 
999  CONTINUE 
ENO  file  9 
END 

SUBROUTINE  RLUT2(FtFMIN*FMAX«NVAR»NFUN«N)iNtXOiOELX) 

PLOT  FIRST  N POINTS  OF  UP  TO  26  FUNCTIONS  FIX) 

F(I.J)  CONTAINS  the  VALUE  FOR  TME  JTM  POINT  OF  THE  ITH  FUNCTION 
FMINlI)  ANO  F*4AX(I)  CONTAIN  THE  HiN  ANO  MAX  ORDINATE  VALUES  FOR 
THE  ITh  FUNCTION. 

NVAR(I)  an  array  of  TITLES  FUH  THE  VARIOUS  FUNCTIONS 
TO  BE  PLOTTED  AOAINST  THE  ABSCISSA 
NFuN  NUMBER  OF  FUNCTIONS  TU  BE  PLOTTED  • DIMENSION  OF 

NvARi  FMIN.  FNAX 

N1  USED  ONLY  IN  F(Nltl)  AS  PASSED  UIMENSION 

N NUMBER  OF  POINTS  IN  A SINOLE  PLOT  FRAME 

AO  first  ABSCISSA  VALUE 

OELX  ABSCISSA  INCREMENT 

DIMENSION  9STEP (26) .F (NI.N) .FMININFUN) tFMAX (NFUN) iVLAST(26) * 

1 VFIoST'(26),HEA0|6)*STCH(26> 

INTEOER  Ch(26) .NVARI  NFuN) .UOT.ASTER.PLUS. BLANK 
INTEOER  C 
integer  a (191) 

DATA  BLANK, DOT. ASTER.PLUS/lH  tlH.tlH*(lH*/ 

DATA  CH(1)  •CH(2)  .CHO)  »CH (A)  .CHiS)  .CH (9)  tCH(7)  .CHtB)  iCH(9)  .CM 1 10) 

2 / |MA  , IHB  • IHC  * IHD  , iHE  • IMF  . liO  • IMH  tlHl  *1HJ  / 
DATA  CHIll) ,CH(12) ,CH(I3|*CH(1A),CH(15) ,CHI16) .CHIlTl .CH(IS) 

2 / IHK  , IHL  t IHH  IMN  « IHO  , IMP  . IHO  . IHR/ 

DATA  CM (19) ,CH(20)«CH(21)tCH(22)«CH(23) |CH(2A) tCH(2S) ,CH(26) 


COPY 

MAIN  179 
MAIN  160 

main  111 

MAIN  112 
MAIN  113 
MAIN  IIA 
MAIN  IBS 
MAIN  1B6 
MAIN  187 
MAIN  IBS 
MAIN  1B9 
MAIN  190 
MAIN  191 
MAIN  192 
MAIN  193 
MAIN  194 
MAIN  19S 
MAIN  196 
MAIN  197 
MAIN  19B 
MAIN  199 
MAIN  200 
MAIN  201 
MAIN  202 
MAIN  203 
MAIN  204 
MAIN  20S 
MAIN  206 
MAIN  207 
MAIN  20B 
MAIN  209 
MAIN  210 
PLOTI  2 
PLOT!  3 
PLOTI  4 
PLOTI  S 
PLOTI  6 
PLOTS  7 
PLOTI  B 
PLOT2  9 
PLOT2  10 
PLOT2  11 
PLOT2  12 
PLOT2  13 
RLOT2  14 
PLOT2  IS 
PLOT2  16 
PLOT2  17 
RLOT2  IB 
HLOT2  19 
PL0T2  20 
RLOT2  21 
PLOT!  22 
PLOT2  23 
RLOT2  24 
PLOT!  25 
PLOT2  26 
PL0T2  27 
PLOT2  20 
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V#.,. 


B!SI  AVAIUBIE  COPY 


I 


I 

I 


Sit 

y 

f 


-d. 

I 

i 

f 


[■  Vi, 


2 / IHS  I IHT  • IHO  * IHV  I 1H«  • IHX  • IHT  , IHZ  / 

c 

ircNruN.tc.e.OR.N.LC.oi  rzturn 
C RRiNT  HCAOINas. 
wRire<««s6) 
kb  rORMAT  »///l 
DO  AO  I«li»4ruN 

30  TCNMaAdSCrMAXCn-RMlNtni 
CXP»». 

IP  (TCNM.rO.O.I  00  TO  2 
C RRINO  TCNM  TO  A vAluC  SCTNICN  1 AND  10 
IPATCNM.LT.l.)  00  TO  1 
3 IPATCNM.IT.IO.)  00  TO  2 

exR»6XP*l^,  ‘ 
riNH«TeNM*.i 
00  TO  3 

1 CXRsCXRs.t 
TENMaTtNM*10. 

IPtTENM.OT.U)  00  TO  2 
00  TO  1 

C SET  UR  VALUE  mCTwEEN  0R(0  LlRESi  RSTtR, 

2 RSTBRsS. 
tP(TENM.ac.S.)RSTCR>10. 
lP(TENH.LT.t*)RSTER>2. 

5 RSTER(I)a<>sTER*EXR**l 

c cumrute  value  op'startino  line*  vpirst. 

riRST«PMlN(I)/RSTER(l> 

IP(PMIN(li,LT,0.»PlRlTaPlRST-i. 

P1RST>AINT(P1RST) 

VPIRST  (llaPlRSTiRSTERdl 
C CHECK  ENO  line  VALUEiVLAST. 

VLAST<I)avPlRST(I)*10.*RSTERil> 

IP(VLAST(n.OT«PMAX(lMaU  TO  « 

C IP  ORARH  IS  TUO  SMALL  TAKE  NEXT  LAROtR  STEP* 

AAsRSTtP 

iPUA.LT.n.iRSTERsS. 

IP<AA«eO.K*)PSTERalO« 

IPCAA.Li.tO.t  00  TO  B 
RSTCR-2. 

EXRalO.«ExP 
00  TO  » 

C COMPUTE  VALUE  OCTaCCN  ROINTS.STE^-. 
k STER(I)aR»!TER(II»,l 
RK«0. 

00  6 P.K«li6 

HCA0<KM "WPIRST ( 1 1 *2,*RK*R8TtP< I) 

6 RKaRK*!. 

AO  RRITE  (6«i*SI  CHdIt  NVAR<n*  <HCAU(KK)tKKal*6l 
AS  P0RMAT<U*A1*3H  a tA10t9XilPEl2*A*S(eXf IPE12*A) ) 

00  50  J>l*tOI 
A(U)aBLANK 

lP(MOO(J«l'i)  *£0.11  A(J)aOOT 
SO  CONTINUE 
' WRlTC(6f9S)  A*A 

SS  PORMAT  (25X*10Ul/l9X*AHTlMe««Xf lOlAl) 

C PLOT  EACH  POINT 
00  too  Jal.N 
BaX0*PLOAT(J-l)«0ELX 
00  TO  KalilOl 


ML0T2  29 
HL0T2  30 
PLOrZ  31 
PL0T2  32 
PL0T2  33 
PL0T2  3A 
PL0T2  35 
PL0T2  36 
PL0T2  37 
PL0T2  3B 
PL0T2  39 
PL0T2  AO 
PL0T2  Al 
PL0T2  A2 
PL0T2  A3 
PL0T2  AA 
PL0T2  A5 
PL0T2  kb 
PL0T2  A7 
HL0T2  AS 
PLOTS  A9 
PL0T2  50 
PL0T2  51 
PL0T2  52 
PL0T2  S3 
PL0T2  SA 
PL0T2  S5 
PLOTS  56 
PLOTS  57 
PLOTS  5S 
PLOTS  59 
PL0T2  60 
PLOTS  61 
PLOTS  62 
PL0T2  63 
PL0T2  6A 
PLOTS  65 
PL0T2  66 
PLOTS  67 
PLOTS  61 
PLOTS  69 
PLOTS  70 
PLOTS  71 
PLOTS  72 
PLOTS  73 
PLOTS  7A 
PLOTS  75 
PLOTS  76 
PLOTS  77 
PLOTS  78 
PLOTS  79 
PLOTS  80 
PLOTS  81 
PLOTS  82 
PLOTS  63 
PLOTS  SA 
PLOTS  85 
PLOTS  06 
PLOTS  87 


I 

1 


•i 

‘1 

1 


V 

• » 

;i 


■S 
• 1 
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la*a*.^ua  krfa; 


o ono  floor*  c»r>  oor>  r>or» 


m AVAIlABlt  COPY 


AiKXatANK 

inMOOIKilO).tQ*l)  A(K)aOUr 
ir(MQO(JiS  A<K»«0UT 

70  CONTINUE 
■ 00  *0  !■  ,NFUN 

L0C«  < (F ; I , JJ -VF IRST  < n > /StC^* 1 » ♦ i .5> 

CaA(L0C' 

A(L0CJ"CH(n 

IF (C.NC. BLANK, AND. C.NC.OOT)  AtLOO-ASTCR 
60  CONTlNUl 

iF(Muo(j*ie).eci.i»oo  fo  os 
WIIIYe<6,SS)  A 
NS  format  (BSXtIOlAlt 
00  TO  100 
OS  MRlTe(0«l4)6«A 
IS  FORMAT  (l2XtlRClt.4..U*10lAl) 
too  CONTlNUl 
RETURN 
CNO 

SUBROUT  I Nl  KUTMCR (NO .T .M ♦ TO  tCFSC ,A,rtC*»FlR8T) 

DIMENSION  vO  (01  iTl  (0)  lYtlO)  iFOlA)  »Fl  (A>  iF2(6)  lEFSKO)  *A(A) 

COHMON/UUT /NPR I NT  «NRLOT lEND 

COMMON  /ACCEL  / XACCLiBOACLiCOACLiHL 

DATA  NAMUMAM2  /IHYlitHVl  / 

NO  ■ number  UF  EQUATIONS*  NO.  UF  COMRONENTS  OF  TO 
r > INHERENUENT  VARIABLE 

H ■ increment  for  MHICM  solution  is  TO  be  RETURNtO  ♦ OR  - 
VO  • THE  VECTOR  OF  OEFINOINT  VARIABLES.  ENTER  ilTM  INITIAL 
VALUES  at  T and  RETURN  *ITH  VALUES  AT  T*M 

ERSE  • RELATIVE  ERROR  CRITERION  FOR  COMRONENTS  OF  VO  ,07  AB8(A> 

A ■ ABSOLUTE  ERROR  SRiTERlON  FOR  COMRONENTS  OF  v6  *LT.  ABS(A) 
NOTE—  ERSk  AND  A MUST  BE  SR|CIFUO  FOR  EACH  COMRONINT  OF  THE  SYSTEM 
HC*  • THC  smallest  BTlR  size  USEO  in  the  ItfTEOjJATlON 
first  SMOULO  re  0 VMkN  KUTMER  IS  ENTERED  FOR  7Hl  FIRST  TIME 
AFTER  that  FIrIt  IS  I IF  KUTMEm  IS  ENTERED  VI TM  THE  SAME  M OR 
IF  IT  IS  ENTERED  MITM  A CHANOEO  M 
IF  first  is  I the  ERROR  SRiTERlA  CANNOT  BE  MEET  AND  TmE  STER  SIZE 
REDUCED  TO  H/ltB« 


plotz  as 

PL0T2  89 
PL0T2  90 
PL0T2  91 
PL0T2  92 
PL0T2  93 
PL0T2  9A 
PL0T2  9S 
PL0T2  96 
PL0T2  97 
PL0T2  98 
PL0T2  99 
PLOT2100 
PLOT210I 
OL0T2102 
MLOT2103 
PLOT210A 
•'LOTZIOS 
PLOT2106 


2 

3 

A 

5 

6 

7 

8 
9 


IF  (FIRS'.) 


10,10,20 

- - - FIRST  rNTHY 


10 


MC  ■ M 
IRLOC  « 
FIRST  • 


1 

1. 


OTHER  ENTRY 


20  LOC  " 0 
MCA  ■ HC 

IF  (MC.NE.I.l  00  fO  30 

BOO  FORMAT*SX?ASMKUTMtR  ENTERED  VITf  ZERO  INTCORATION  INTERVAL 
FIRST  > 2. 

RETURN  ^ , 

......  .....  S CALLS  rO  DAUX 

30  CALL  OAUX(T,VO,FO) 

■ 1F(NPR1NT,!O.8)fR17E(6,A00)V0»T,F0 
ABO  FORMAT  (6(2X,F10.A)  *AHT1mE*2x,F  lO.o^) 
IF(NPRINT,FO,8|*RITE(6*AOO)KC 
39  DO  AO  IM.NO 


KUTMER 
KuTMER 
KUTMER 
KUTMER 
KUTMER 
KUTMCR 
KUTMER 
KUTMER 
KUTMERIO 
KUTMERII 
KUTMER12 
KUTMER13 
KUTMERIA 
KUTMERIS 
KUTMER 16 
KUTMERIT 
KUTMERIO 
KUTMCR19 
IKUTMCR20 
KUTMER21 
KUTMER22 
KUTMER23 
kutmehza 
kutmerzs 

KUTMER26 

•'UTMERZT 

i^UTMERZa 

KUTMER29 

KUTMER30 

KUTMCR31 

KUTHeR32 

KUTMCR33 

KUTMCH3A 

KUTMCH36 

KUTMCR36 

KUTMCR37 

KUTMER3S 

«.UTMCR39 

KUTMCRAO 

KUTMERAl 
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BEST  AVAILABLE  COPY 


C 


c 


c 


c 


c 


c 

c 


c 


c 


%o  riLU  ■ vO(i)*(pic/3*)*ro(i> 

ir(NPRlNT,CO.S)riRlTC<6f«00)vt.T 


KUTMER42 

KUTMEB03 


CALL  0AUX(T*HC/3««Yl«Pn 
If (NPRlNT.fO.S)riRIte(6*400)ri*1 
00  SO  I>lt'AO 

SO  ykil  ■ YO(n«(Hc/6.)»fo<i)*(HC/b.i*ri(i» 

If  <NPRlNT,rO.S)MRne(OfAOO)Vl*T 

CALL  0AUX<T*HC/3»fYlifit 
If (NPRINT.'O.S>HRtTe<O*40»)ftit 
00  00  1*1 fNO 

00  Ti(i>  ■ YOtn*<He/o,»»fo<J>*.3Ts«HC»fni» 

If  (NPRlNT.eo.S triRITC <0*400)  YUT 

CALL  0AUX(T*HC/2.fYlifl) 

If <NPRINT.'O.S)HHlTe<0*400)f2tT 
00  TO  I*1«N0 

70  YlLl)  ■ yO(I)«<HC/l.)«fO<t)-I,S»MC»f|<l)**.»HC»f2<I) 

If  <P<PRINT,'O.S)PRlTC(0»40f>YliT 

CALL  OAOX(T*HC*Yl«fl) 

If <NPRlNT.'O.S)«RITe(O*40O>fttT 
00  00  I>t«NO 

00  Y2(I)  ■ Y0(I)*HC/O,«f0<I)«<ti/3.)*»4C*f2(I)*<HC/O.)»fl(I> 

If (NPRINT,eO.S)MRITe (0*400) Y2*T 

INC  ■ 0 

..........  CHECK  ERROR  CRITERIA 

00  110  l«l.NP 

222  • AOS(VlLl))-A<I) 

If  (222)  OSiOTiOT 

..........  ariolutc  error 

01  ERROR  ■ A0X(*2t<Yl(I)-Vt<l))) 

If  (CRROR.td))  100*100*00 

..........  relative  error 

07  error  ■ Aes(.2*.2»Y2(I)/Yi(I)) 

If (eRROR.CBSC(l))  100*100*00 

...........  SINCE  error  .OT*  error  CRITERIA  CHECK  IP  HC. 

..........  IP  YES  THEN  HALVE  INTERVAL*  OTHERWISE  STOP* 

00  X ■ I20.*AhS(HC).ABS(H) 

If(X)  Ol*OS*OS 

..........  error  TOO  LAROE 

01  WRITe(«*V2iI*T*ERR0RiHC 

02  PORRAT(/ISh  pur  equation  NO.  12*27***  THE  RELATIVE  ERROR  AT  T 
t ElS.Oi  AH  IS  tElS.OtlSH  STEP  SIZE  ■ *£15.0) 

PIRST  ■ 2, 

RETURN 


KUTMCR4A 
KUTMCR45 
KUTMERAO 
KUTMER47 
KUTHER40 
KUTHER44 
KUTMERSO 
KUTMERSl 
KUTNCRS2 
KUTHERSS 
KUTHCRS4 
RUTMERSS 
KUTMERSO 
KUTMERS7 
KUTMERSO 
KUTMERSO 
KUTMEROO 
KUTMEROl 
KUTMCR02 
KUTMffiOS 
KUTMEROO 
KUTMEROS 
KUTMEROO 
KUTMEH07 
KUTMEROO 
KUTMEROO 
KUTMCR70 
KUTMERTI 
KUTMER72 
KUTMCfl73 
KUTMERTA 
KUTMCR7S 
KUTMERTO 
KUTMER77 
KUTMERTO 
0T.H/KUTMER70 
KUTMERSO 
KUTMERSl 
KUTMERS2 
KUTMERS3 
KUTMERSA 
■ * KUTMERS5 
KUTMERSO 
KUTMEMST 
KUTMERSO 


..........  *1ALVC  INTERVAL 

OS  HC  ■ HC/2. 

IPLOC  ■ 2*IPLUC 
LOC  » 2»L0C 

hcx  ■ hc 

WRITE(2f7|K)T*I*CRR0R<HC 

TIO  PORMAT(/BH  TIME  ■ <PI0.3*SX«20HHALVE  INTERVAL.  EQUATION  *13* 
'•13H  HAS  error  ■ |EI0.0*0X*I7H  STEP  SUE  NOW  ■ «E1S.S) 
WRITC(2«72K)  NAM2.(r2<J)«JRltN0) 

WRITEIZtTZK)  NAMI*(Yl(J)*J«).ND) 

720  PORMATI  2X.A2  / 3<10E13.S/)) 

00  TO  30 


KUTMERSO 

KUTMEROO 

KUTMEROl 

KUTMER02 

KUTMEROO 

KUTMEROA 

KUTMEROS 

KUTMEROO 

KUTMER07 

KUTMEROO 

KUTMEROO 

KUTMEIOO 


56 


BEST  AVAIUBIE  COPY 


C • * 

c 


.........  test  ir  INTCHVAt  LENOTH  CAN  BE  DOUBLED 

> ir  (ERftOR«6A..E^SE(I)l  UO«tlO»lDl 
i INC  • I 

> CONTINUE 

.........  UPDATE  t AND  SOLUTION 

I T • T*MC 
DO  at  I>\,ND 
i YO(I)  • Vt(I) 

.........  MT  SOLUTION  IN  NEXT  INTERVAL 

LOC  X LOC*l 

If  (LUC'IPLOC)  iio<tio«aio 

> If (INOaiO.IDO.lIO 

If  (LOC-(LOC/^>*2)  >10«140*2IO 
} If (IPLOC-I)II^2lO«m 
.........  double  interval  lenoth 

) HC  • 2,*HC 
LOC  ■ LUC  /t 
IPLOC  ■ IPLOC/f 
) If <IPLOC-LOC) 

) BbACL  ■ fO(2)-MACCL*fO(3l 
COACL  • f0(2l 

return 

END 

END 

SUBROUTINE  DAUX(T:mE*X«RHS> 

TIME  TIME  AT  WHICH  SYSTEM  IS  TO  BE  EVALUATED 

X STATE  VECTOR 

RHS  the  RIOHT  hand  SlOE  Of  THE  EQUATION  S • f A 

REAL  RAM 

REAL  IA«IT.M«X*MA*MASS«NC(»«NL*NtMMAX 
INTEOER  CNO«PTtME 

dimension  < (6> «RHS(«) if 0»)T»A(3»3) t index 0»3»« 
t R<120l  iVdaOTiOtUO) 

common  /ship/  MASSiClNTiQAiCE*CE2iCC3iDMUiC0MUiE20MUiE30MUiBfiBMM 

• NLifLiIA.l(130) 

COMMON  /CONST/  NC6fEC0iPIi0PR«RPDiORAVTYiRHOiKiNUM|MA(l20)iC0iTAi 

• B<120)  iSETAiHWdlOTiTIiUNAUiWtXDiTiXPiM.lTi 

• OrLTASiYXtESTdaOtiCiROfHARiMMAXd.  0>iTESTd20)  i 

• Nd20)iPHALf 
COMMON  /IN/  nM(|20).Bld20>,VELIN 
COMMON/UUT /NPR I N iNPLOT  lE  1^0 
COMMON  /SEAnAVE/  STARTiRISE.RAMP 

COMMON  /vAvE/  PiPT(120)i2HA,Z«MA*LHASdlWMA,ZBEMAi22MMAiE2MAZi 

• 2B0nT<120) 

ramp  ■ HMPiTlMCiSTARTiRlSl) 

PIM  ■ Pl/a, 

CT  ■ C*TIME 
CX6  ■ COS d <6)1 
SX6  ■ SlNtx(6) ) 

»*t*SCT  VALUES  Of  MA  AND  B 
DO  TB  l»l,NL'M 

PTd)  • (X(A)*Ed)*CXeiNd»»S»R*CT»*R 
Rd)  • «0*COS<PTd)l«RAMP 

!••••••  COMPUTE  HW  SUBMEHOENOk  Of  A POINT  AND  R THE  WAVE 

NW(I)  IS  IN  THE  7mO  COOPUINATE  SYSTEM 


KUTMElOl 
RUTMEI02 
KUTME103 
KUTMEIOA 
KUTMEIOS 
RUTME106 
RUTME107 
KUTMEIOB 
RUTME109 
KUTMEllO 
t^UTMEUl 
KUTMElia 
KUTME1I3 
KUTMEIIA 
KUTMEIIS 
KUTME116 
RUTMEU7 
MUTMEIIB 
RUTME119 
KlTMElaO 
KUTME121 
KuTMEiaa 
KuTMEiaa 
MUTME124 
UAUX  2 
OAUX  3 
OAUX  4 
OAUX  S 
OAUX  6 
OAUX  r 
OAUX  S 
OAUX  0 
OAUX  10 
OAUX  II 
OAUX  12 
OAUX  13 
lOAUX  14 
OAUX  IS 
OAUX  lA 
UAUX  17 
OAUX  IS 
OAUX  19 
OAUX  20 
UAUX  21 
UAUX  22 
OAUX  23 
OAUX  24 
UAUX  2S 
UAUX  26 
OAUX  27 
UAUX  2S 
UAUX  29 
UAUX  30 
OAUX  31 
OAUX  32 
OAUX  33 
UAUX  34 
UAUX  39 
UAUX  36 


BEST  AVAIIABIE  COPY 


MM(ll  R X(M)-E<I)*SX6*N(n*CX6-R(l) 

UAUX 

37 

iFEhwtn.ar.O)  00  to  «s 

OAUX 

3a 

Craft  is  not  sushcrolo 

OAUX 

39 

MA<h  * 0. 

OAUX 

40 

eitixo. 

OAUX 

41 

0(1)  ■ 0. 

OAUX 

42 

00  TO  T4 

UAUX 

43 

6S 

V(l)  ■ -R0»K*SIN(RTI1))*RAMR 

OAUX 

44 

0(1)  • H«(|)/(CX6-V(I)*|XA) 

OAUX 

4S 

D(t)  is  IN  THC  lOOV  AXIS  SYSTCM  AND  IS  THE  SUSMERbENCC 

OAUX 

4S 

1F<0(1).0E.TCST(1))  00  TU  TO 

UAUX 

47 

Craft  is  partly  susmehoeu 

OAUX 

40 

H<I)  ■ 0<I)«(1,/TA)*P1H 

OAUX 

49 

11(1)  ■ 0(II*(1,/TA)*P1H 

UAUX 

50 

NA(I)  ■ KAR«PMALF*R(I)*B(1) 

UAUX 

SI 

00  TO  7S 

OAUX 

S2 

CHINE  IS  IMMERSED 

OAUX 

53 

01  ARRAY  IS  USED  FOR  THE  INTEORALS  OVER  THE  PORTION 

OAUX 

54 

OF  THE  HULL  FOR  MHICH  THE  CHINE  II  NOT  IMMERSED 

OAUX 

SS 

MA(1)"MMAX(II 

OAUX 

56 

li(I)>IIM(I) 

OAUX 

57 

U1(I)>0* 

OAUX 

58 

71 

continue 

UAUX 

59 

IF (NPR1NT.lt. 41  00  TO  IS 

OAUX 

60 

SRlTE(Sf74)TlHE 

OAUX 

61 

7A* 

FORMAT)"  TIME  ■ 'SF10.4) 

OAUX 

62 

NR1TC(S«7A)  U(l)«IalfA) 

OAUX 

63 

MRlTEtStT?)  (R(l) iIrIiNUM) 

OAUX 

64 

WRITE  (A«7R)  (HS(l)«IalfNUM) 

OAUX 

6S 

WRITE  (StTR)  t B(l) •I■I«NUM) 

UAUX 

66 

WRITE  (V  (1)  tlaltNUM) 

OAUX 

67 

WRITE  (SiStI  (0(l)«l*l«NUM) 

OAUX 

61 

NRITE(S«St) (MA(I»ilaliNUM) 

OAUX 

69 

7» 

FORMAT)"  X(I>  "lAdKiElltOI) 

OAUX 

70 

77 

FORMAT  ("  R(l)"f lOFlOiA) 

OAUX 

71 

7« 

FORMAT  ("  HW)l)''*l«F\0.4) 

OAUX 

72 

7W 

FORMAT  ("  R(I)»«10F10.4) 

OAUX 

73 

00 

FORMAT  ("  v(I)"il0F10.4) 

OAUX 

74 

•I 

format  ("  O(l)"<l0F10*4) 

OAUX 

75 

■e 

FORMAT)"  MA)I)  "(lOFlO^A) 

OAUX 

76 

AS 

continue 

UAUX 

UAUX 

77 

70 

• • 

••••••  COMPUTES  NL  ANO  FL  ANU  THE  ASSOCIATED  INTEROALS 

OAUX 

79 

CALL  FuNCT)XI 

OAUX 

OAUX 

00 

61 

if)NPrint.lt.4)00  to  it 

UAUX 

82 

WRITE  )6i IS)  TX.FLiORAO.TJfW.NL.XOtT.XP 

UAUX 

03 

IS  FORMATE"  ".lOeiZ.S) 

OAUX 

04 

IT 

CONTINUE 

UAUX 

as 

• • 

a • • • • COMPUTE  THE  F VECTOR 

UAUX 

06 

F|ltl)  ■ .TX*FL*SX6«0RA0«CX« 

OAUX 

07 

Fd.DaO.O 

OAUX 

60 

F(lfl)  ■ TZ«FL«CX6t0RA0*SX««W 

OAUX 

09 

FI3f l)aNL-nRAO*XO*T«XP 

OAUX 

90 

IF)NPRINT,LT.3iao  TO  IS 

OAUX 

91 

WRlTOSilO)  )F)I.l)«Iali3) 

OAUX 

92 

IS 

continue 

OAUX 

93 

t • 

« • • • • COMPUTE  THE  A MATRIX 

OAUX 

94 

A)ltl)  ■ M*MA5S*SX6*SX6 

OAUX 

9S 

58 


BEST  AVAIUBLE  COPY 


A(l<2)  ■ MASS»SA6«CX6 
A(l,3)  ■ -()A«iX6 
A<1<2>  ■ 0. 

A(l<3)  ■ 0. 

A(2«1)«A(1,2) 

A(2t2)  ■ M*HASS*CA6»CX6 
A(2t3)  ■ -OA^CXft 
A(3tll<iA(t.3l 
A(3«2)aA(2«3) 

AO«3)alT*IA 

ir(NPRINT,LT.3)00  TO  23 
MRlTe(6il2)  (A<l«U«la)t3l 
MRlTe(6*l3)  (A(I«2)tlaU3) 
HRire(6iU)(A(l«3I.I>l,3) 

invert  the  a matrix 

23  CAUL  MATlN^IA'OOtXtltltOETERMtlUtlNDEX) 

ir  <lt).EQ.2)«RITE<6«2«) 

26  FORMAT  ("  MATRIX  1$  SlNOUtAH  ••) 

C»««««*A  ON  RETUON  mill  CONTAIN  TmE  INVERSE  MATRIX 
C IO>2  MATRIX  Is  SINGULAR 

C ■!  INVERSE  MAS  FOUND 


!••••••  COMPUTE  THE  right  MANl)  SIDE 

RH8(l)  ■ F(l,l) 

RHS(2)  ■ F(2«n 
RHS<3)  ■ FI3<1) 

Rh3(1)  " 0.0 
RHS(4)  ■ X(l) 

RHSO)  ■ X(2) 

RHS(6)  ■ X(3> 

10  FORMATS  F(I,n  •S3(2X<C12.M)| 

12  FORMATC  AUtl)  ••t3<2X*Et2.An 

13  FORMATS  A<Ii2l  ••.3<2XfEI2.«> » 

lA  Formats  Ad.s)  ••.3(2x«ei2.ah 

30  IF (NPRINT,LT.2I  00  TO  AO 

HRlTE(6tl2t  (A(t«l) ilMltS) 

•RtTE<6d3)  (A(|,2)  tl*lt3> 

MRITE(60AI  (A<I»3I  ilMlOl 
HRITE<605)  (RHS(II  tlaltA) 

33  FoRMATC*  RHS(I)  ••«6(2XiE12.6I  ) 

AO  CONTINUE 
RETURN 
ENO 

SUHROUTINE  FUNCT(X) 

REAL  KAH 

REAL  lAtlAA.IPARTiKtKPIiMA.MASSfNU.NCO.lT.MtMMAXtN 
INTEGER  END 

DIMENSION  IPART(I20)«C1 U20| ,C2 (120) , 

t 01(120) i02(120l >03(120)*UA()20) iDS ( (20) *06 ( 120) t 

• OPART(120)iZl(U0)  t^2(UO>.23()20)  .ZAdEO)  iZ5d20)  • 

» 26(120)127(120) 

• .X(6)  iVHAAdEO) 


DAUX  96 
OAUX  97 
OAUX  90 
OAUX  99 
OAUX  100 
OAUX  101 
OAUX  102 
OAUX  103 
OAUX  lOA 
OAUX  103 
OAUX  106 
OAUX  107 
OAUX  106 
OAUX  109 
OAUX  110 
OAUX  111 
OAUX  112 
OAUX  113 
OAUX  llA 
OAUX  113 
OAUX  116 
OAUX  117 
OAUX  116 
OAUX  119 
OAUX  120 
OAUX  121 
OAUX  122 
OAUX  123 
OAUX  126 
UAUX  123 
OAUX  126 
OAUX  127 
OAUX  126 
OAUX  129 
OAUX  130 
OAUX  131 
OAUX  132 
OAUX  133 
OAUX  136 
OAUX  135 
OAUX  136 
DAUX  137 
UAUX  136 
FUNCT  2 
FUNCT  3 
FUNCT  6 
FUNCT  3 
FUNCT  6 
FuNCT  7 
FUNCT  8 
FUNCT  9 
FUNCT  10 


FUNCT  11 

COMMON  /SHIP/  MA9StClNT«OAtCE<ce2tCE3.DMUtEOMU«E2OMUiE30MUiBF»eMM«FUNCT  12 

• NL.FL,IA«Ed20)  FUNCT  13 

COMMON  /CONST/  NCO.EC0tPIiOPR>(<.>O«OHAVTT»RHU,K«NUMtMA(120)  tCOtTAi  FUNCT  U 

• 0(120)  tBETAiHMdZO)  tTZiOMAUiMiXOtTtXP.MdT,  FuNCT  13 

• 0rLTAs«TX,ESTd20),CtRU.KXR,HMAXd  0)  «TEST  ( 120)  • FUNCT  16 

• N(1Z0).PMALF  FuNCT  17 


BESr  AVAIIME  COPY 


COHMON  /IN/  SM(120) «Bl<ll9>«VtLlN 

PUNCT 

10 

CONHON/UIJT  /NPRI  NT  »NPCOT  tCNO 

PUNCT 

19 

COMMON  /wave/  R(120|.PT(l20).2MA,2wMA»CMAS»Z2WMAtZwEMA,22«HAi 

ezmazpunct 

20 

« 

•ZW00T<120) 

PUNCT 

21 

COMMON  /inter/  IItKTT(10)«01PF<10) 

PUNCT 

22 

COMMON  /SEAWAVE/  START *R1SE«RAMP 

PUNCT 

23 

COMMON  /TEAT/  VMA 

PUNCT 

26 

• 

• •••••••  initialize  inteoral  sums 

PUNCT 

2S 

MASS  ■ O.O 

PUNCT 

26 

OA  ■ O.O 

PUNCT 

27 

lA  ■ O«0 

PUNCT 

20 

CC  ■ 6t0 

PUNCT 

29 

CC2  • o;o 

PUNCT 

30 

OMU  ■ t.O 

PUNCT 

31 

EOMU*0.0 

PUNCT 

32 

E20MU'*  0.0 

PUNCT 

33 

ESOMU  *0.0 

PUNCT 

34 

tfP  ■ 0.0 

PUNCT 

3S 

BMH  a'O^O 

PUNCT 

36 

2MA  ■ 0.0 

PUNCT 

37 

ZWMA  >•  0.0 

PUNCT 

36 

CMAS  ■ 9.6 

PUNCT 

39 

ZZWHA  a'O.O 

punct 

40 

ZWEMA  a 0.0 

PUNCT 

41 

22WMA  a 0.0 

Punct 

42 

ESMAZ  a 0.0 

PUNCT 

43 

VPART  a X(U»SIN(X<6n*X(2)aCUS<A(«>> 

Punct 

44 

SX«  a S1N<X(6)) 

PUNCT 

4S 

CXA  a COS(X<«>> 

PUNCT 

46 

WO  a KaC 

punct 

47 

• 

• a a • • • SET  UP  THE  PUNCTIUNI  TOR  THE  INTEORALS  (PAOE  A 

OP  NOPUNCT 

4S 

00  90  laUNUH' 

PUNCT 

49 

IPART(l)aE(I)aE<I)aMA(I> 

Punct 

SO 

OPARTCnaCdiaMAd) 

Punct 

SI 

ZWOOTd*  a -RoawO*SIN(PTd)»aMAMP 

PUNCT 

92 

U a X<liaCiA-X(2|aAX6a2iOUTd)aSX9 

PUNCT 

93 

VEL  a vHARt.XOiaEdl-ZwOUTdiacxp 

PUNCT 

S4 

Zld)  a MAdtaZWOOTd) 

PUNCT 

99 

Z2d)  a -MACDacOSlPTdliaRAMP 

PUNCT 

96 

Z3d)  a EdiaZZd) 

Punct 

S7 

ZAd>  a Ediazid) 

PUNCT 

SO 

ZSd)  a uaZ2d) 

PUNCT 

59 

26 d)  a C(n*Z9(n 

Punct 

60 

ZTdl  a MAdiavELau 

Punct 

61 

ir  (VEL.LE.O.)  00  TO  60 

Punct 

62 

ir  (Rid) .le.o.oi  00  to  so 

Punct 

63 

OROT  a Za00T(l)a(Xd).C»X(3)a(Nd»aCX6-ed>*SX6)»/C 

Punct 

64 

Old*  a vELaSl  dta(X(2)«X(3)a(CX6aEd>*SX6aNd)  1 -OROT) 

Punct 

6S 

00  TO  SI 

punct 

66 

5U 

01d>  a 0. 

PUNCT 

67 

Si 

CONTINUE 

PUNCT 

60 

02d>  a EdiaOld) 

PUNCT 

69 

Cld)  a vCLaVELaHd) 

PUNCT 

70 

C2d>  a E(naCld) 

PUNCT 

71 

00  TO  61 

Punct 

72 

SO 

Old)  a 0, 

PUNCT 

73 

02(1)  a 0. 

Punct 

74 

Cld)  a 0. 

punct 

75 

C2d)  a 0. 

punct 

76 

BEST  AVAIWBLE  COPY 


61  CCNTlNlje 

03 (I>  • Z2(I)*VEL 
06 <11  • C(I)«03(|) 

PIM  • PI/2. 

.05(1)  ■ B(|)«(Htt(I)-B(I)«IA/2.) 

66  06(11  ■ 05(I>H(n«.S 

90  COMTINUE 
RHOaar.H0*GUAVTV 

C ••#••••  sCT  UP  TMC  FUNCTIONS  FoN  THE  INTEOPALS 
PIM  m Pl/2. 

KPI  > KAR«PI 

C evaluate  INTEOPALS  USING  TPAP  METHOU 
I • I 
INDEX  > I 

91  CALL  TPAP(HA(IN0CXI,0IFF(1),KTT(1).TMASS) 

CALL  TRAP (OPART ( iNOEXl f OIFF ( 1 1 *KTT ( 1 ) iOaI ) 

CALL  TRAP (Cl ( INOEX) <OIFF  <1 > *kTT ( I > fCEA) 

CALL  TRAP (C2( INDEX) tOIFF < I ) *KTT ( 1 ) f CC2A) 

CALL  TRAP ( IPART (INDEX ) tOIFF ( I ) (KTT ( 1 > • 1 AA ) 

CALL  TRAP  (01  ( iNOEX)  lOIFF  ( 1 ) ,<(TT  ( 1 > «UMuA) 

CALL  TRAP (02 ( INDEX) tOIFF ( 1 ) .KTT ( 1 ) •EDHUA) 

CALL  TRAP (03 (INDEX) tOlFF (I ) .KTT ( 1 ) tiZDHuA) 

CALL  TRAP (06 (INDEX) lOlFF ( 1 > .KTT ( I ) tC3DMuA) 

CALL  TRAP (nS( INDEX) tOlFF (I ) |KTT ( 1 ) tttFA) 

CALL  TRAP (06 ( INDEX) lOIFF ( I ) .KTT ( 1 ) •BHHA) 

CALL  TRAP (Zl (INOEX) «OIFF( I ),KTT(1)»ZHAA) 

Call  trap (72( index) .oifF(1).ktt(1) «2«naa) 

CALL  TRAP (73 (INDEX) tOIFF (1) ,KTT(1) »EMASA) 

CALL  TRAP(7«(IN0CX) ,01FF ( 1 ) .KTT ( I ) «22«HAA) 

CALL  TRAP ( 75 ( INDEX) »01FF ( I ) .KTT ( 1 ) tZwEMAA) 

CALL  TRAP(76(IN0CX) iOIFF ( 1 ) iKTT ( 1 ) tZEpMAA) 

CALL  TRAP (77 (INOEX) tOlFF ( 1 ) ,KTT ( 1 ) tEZHAZA) 

C 

93  CONTINUE 

HASS  • •’'ASX  * THASS 
QA  • QA  * OAl 
lA  • lA  • lAA 
ce  • CE  * CEA 
CEZ  ■ CE2  * CE2A 
OMU  ■ OH')  • ONuA 
EDMU  • EOMii  * EDHUA 
E20HU  * C20MU  * E20HUA 
E3DHU  ■ E30MU  * E3DNUA 

BF  ■ BF  4 oHUG*eFA 
BHM  ■ BHM  4 RHOG*BNHA 
ZMA  ■ 2MA*ZHAA 
ZVMA  « Z^MXfZvMAA 
ENAS  ■ EHA<-,4EHASA 
ZZHHA  ■ ZZxHA427wMAA 
ZmEMA  ■ ZWEHA42NEMAA 
ZZmHA  3 Z2wNA422hHAA 
C2MAZ  ■ E2HAZ4E2HA2A 

96  CONTINUE 

IF  ( I.OE.IDGO  to  92 
INDEX  • INnEX»KTT(l)-l 
I ■ 1*1 
GO  TO  91 

92  CONTINUE 


FUNCT  77 


FUNCT  78 
FUNCT  79 
FUNCT  80 
FUNCT  81 
FUNCT  82 
FUNCT  83 
FUNCT  86 

(PAGE  5 UF  NOTES) FUNCT  85 
FUNCT  86 
FUNCT  87 
FuNCT  88 
FUNCT  89 
FUNCT  90 
FUNCT  9i 
FUNCT  92 
FUNCT  93 
FUNCT  96 
FuNCT  95 
FUNCT  96 
FuNCT  97 
FuNCT  98 
FuNCT  99 
FuNCTlOO 
FUNCTlOl 
FUNCT102 
FUNCT103 


FIJNCT106 
FUNCT105 
FUNCT106 
FUNCT107 
FUNCT108 
FUNCT109 
FuNCTllO 
FUNCTlIl 
FUNCTU2 
FUNCT113 
FUNCT116 
FUNCTUS 
FUNCT116 
FuNCTlIT 
FUNCTllO 
FUNCT119 
FUNCT120 
FUNCT 121 
FUNCT122 
FUNCT123 
FUNCT126 
FUNCT125 
FUNCT126 
FUNCT127 
FUNCT128 
FUNCT129 
FUNCT130 
FUNCT131 
FUNCT132 
FUNCTID3 
FUNCT136 
FUNCT135 


61 


BtSI  AVAIWBlt  COPY 


C 

c 

c 


• • • * • CALL  COMPUV  TO  flNO  TmL  vALU£  OF  nL  ANO  FL  USING 
THC  VALUES  or  THE  AROvE  INTEGRALS 
CALL  COHPUT<A) 


ir<NR«INT,LT.3)  GO  TO  111 
iriNPRINT.EO.S)  00  TO  lOS 
ir(NPRtNT«CO.A»00  TO  lOO 
WRlTClOfOT)  (IRARTdl  i!al«NuM) 

WRlTe<6«9S)  (ORARTin «I«1«NUM> 

«RlTe(6f9«)  (Cin>«|aliNUN) 

MRlTCIOilOO)  (C2<niI>l«NUM) 

MRITEIOilOl)  <C3(t)»I>liNUM) 

RRire(6il«2)  <01<niI>l«NUM| 

WRireiOtlO))  402<I)«I<ltNUM) 

NRiretOilftA)  (03(1) dal tNUH) 

«RITe(6«l09)  (OAdlilaUNUM) 

MRlTe(0«ie6)  (05(I)iIal»NUM) 

WR1TE(««112)  (00(1) flat tNUM) 

HRlTClOtllS)  (ZKDflaliNUH) 

WRITE(6fIl4) (22(I)fIal,NuH) 

MRlreiOillS) (Z3(I)«|al,NUH) 

HRlTEIOfllO) (ZA(l)fIal,NUN) 

WRlTe(OfllA) (ZS(I)fIal«NUH) 

MRirE(6fll9) (Z6(I)*Ial«NUM) 

MRlTE(Ofl20) (Z7(l)«lal,NUH) 
MRlTC(6itOT)KPt*RHOOfPlH 
190  NRlTe(Ofl09)  MASSfCINTfOAfCCfCetfCCS 
liRlTe(0«l2l)lA 
121  FORMAT (•  lA  atClO.A) 

HRm(6iii0)OMUfeDMu«e2OMu«e3aMUfiir,aHM 
NR ! TE ( 0 • 1 1 7 ) ZM A « ZRM A f CHAS  f Z ZRMA , ZmCMA • ZZnMA  f CfHAZ 
aaaaaaaaaa  FORMATS  •••••###•#• 


90  FORMAT  C 

CFART(I)"i 

97  FORMAT (" 

tPART(l)"i 

90  FORMATS 

(JPA9T(I)»( 

99  FORMAT C 

Cl 

ISO  FORMAT  ('• 

C2 

M 

191  format (“ 

C3 

H 

192  FORMAT  (•• 

nl 

II 

193  format (" 

02 

M 

i9a  format C 

03 

II 

195  FORMAT)" 

OA 

II 

190  FORMAT (" 

05 

II 

112  FORMAT)" 

00 

II 

19F  FORMAT)" 

RPHI 

**|C10 

199  FORMAT)" 

MASS 

•"CC2  — _ 

110  format ("  Onu  OfElO. 
•eto.4«"  bf  'sEio.Ai 

113  F0RNAT(AM  zl  i|0(2X 
no  F0RMAT(4h  22  il6(2X 
US  FORMAT  (4m  Z3  tlO(2X 
110  F0RMAT(4H  24  flO(2X 
110  FORMAT (AH  2S  flO(2X 

119  FORMAT (AH  20  ilOlZX 

120  FORMAT (AH  27  flO(2K 


10(2Xfei0.A)) 

io(2xfeio.A)) 

10(2XfC10«A)) 

lOfZXfClOfA)) 

10(2XfE10*A)> 

10(2XiC10.ah 

iO(2x.eio*AH 

10(2XtE10.A)) 

10(2XfE10fA)) 

10(2XtE10.A)) 

lOdXiElO.A)) 

l0(2XiE10.A)) 

>4t»RHU0  "fClOfAf"  PHIH  'SElOfAt 
At"  CINT  'SCtU.A*'*  OA  »fE10.Ai*'  CE  "ftlO.Ai 
••tClO.A) 

At"  EOMU  "tClO.At"  E20MU  "tElO.At"  C3DMU  ••• 
•’  BMM  'seio.Ai 
ElO.AI) 

ElO.A)) 

ElO.A)) 

ElO.A)) 

ElO.A)) 

ElO.A)) 


ilZ  FORMAT (SH  2MA  .EIO.A.OH  2NMA  tElO.A.OH  EMAS  tElO.A, 

• 7H  ZZNMA  ,C10.A,7H  ZWEMA  ,E10.A.7h  ZZPMA  .ElO.A, 

. 7H  E2MA2  .ElO.A) 


FUNCTI 36 
FUNCTI37 
FUNCTUO 
FUNCT139 
FUNCTIAO 
FUNCTIAI 
FUNCT1A2 
FUNCT1A3 
FuNCTlAA 
FUMCTIAS 
FUNCTI AO 
rUNCTlAT 
FUNCTIAO 
FUNCTIAO 
FuNCTlSO 
FUNCT151 
FUNCTI 52 
FUNCTI 53 
FUNCTlSA 
FUNCTISS 
FUNCTISO 
FUNCT197 
FUNCTISS 
FUNCTISO 
FUNCTI 00 
FUNCTlOl 
FUNCTI 02 
FUNCT103 
FUNCTIOA 

functios 
FuNCTlOO 
FUNCT107 
FUNCTIOS 
FUNCTlOO 
FUNCT170 
FUNCTI 7 1 
FUNCTI 72 
FUNCTI 73 
FUNCTI 7A 
FUNCTUS 
FUNCTI 70 
ruNCT177 
FUNCT17S 
FUNCTI 79 
FUNCTISO 
FUNCTiei 
FUNCTI 02 
FUNCT1B3 
FUNCTlSA 
FUNCTISS 
FUNCTISO 
FUNCTIST 
FUNCTISS 
FUNCT189 
FUNCTISO 
FUNCTIOI 
FUNCT192 
FUNCT193 
FUNCT19A 
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BEST  AVAILABLE  COPY 


C 


c 

c 


c 


c 


FUNCT195 
FUNCT196 
fUNCTl97 
COMPUT  2 
COMPUT  3 
COMPUT 
COMPUT  S 
COMPUT  6 
COMPUT  7 

COMMON  /SHIP/  M*SS<CtNT(OAtCCtCe2»ce3tOMli,CaMU»C2UMiJfe3UHU»BF«BHM,COMPUT  B 
• NL*rL»IA«e<t20)  COMPUT  9 

COMMON  /'CONST/  NCa,CCa«Pt»0PRtRPU*GRAVTY»RHO,K,NUM,MA(120) .CO,TA,  COMPUTIO 


B(120) tBETA»Htt(120) «TZ«URAO,H»XO«T 

»XP.M,lTt 

CONPUTIl 

1 OELTAs«TX«CSTtl20l «CiRO»AAR.MHAX < 1 

0|»TEST(120)t 

COMPUT 12 

N(i20)»PHAL^ 

COMPUT13 

COMMON/OUT /NPrt I NT  *NPLOT •END 

CONPUTU 

COMMON  /TEPMS/  TI tT2.T3tTA«T5tT6iT7.TB 

COMPUT 15 

COMMON  /PAVE/  R<I20ltPT(t20)*/MA«ZMMA«EMASt 

ZZWMAiZwEMA.ZZWMA. 

COMPUTIO 

E2MA2.Z«00T(120) 

COMPUT 17 

COMMON  /TC‘iT/  VMA 

COMPUT IB 
C0MPUTI9 
COMPUT20 

CXB  ■ COS (X (61) 

COMPUT21 

SX6  ■ SIN<X(6)) 

C0MPUT22 

DO  ■ K*C 

COMPUT23 

PIN  - PI/2.0 

C0MPUT20 

KPl  ■ KARPPI 

COMPUTES 

CONSl  ■ RO*HO*DO«CX6 

C0MPUT20 

C0NS2  ■ (KOt«RH0*PlH/TA)/CX6 

COMPUT27 

CONS3  ■ A0«W0«K*CX6»SX6 

COMPUT2S 

CONSA  ■ R0*W0*K«CX6*CX6 

COMPUTES 

TEPMl  ■ X<))*CX6 

COMPUT30 

TERM!  ■ X(PI«SX6 

COMPUT31 

UVNUM  • (X(1)*CX6«(X<2)»ZH00T(NUM) )*SX6)* 

COMPUT3E 

. (X(|)*SX6-X(3)«£(NUM>HX<2>>2ND0T(NUM))*CX6) 

COMPUT33 

C0MPUT30 

ZMA  ■ ZMA*X(3)«SX6 

COMPUT35 

22WMA  ■ 2ZWMA*X<3)*SX6 

COMPUT36 

ZWMA  ■ ZPMA»CONSl 

COMPUT37 

ENAS  > EMAS»CON5l 

COMPUT30 

OMU  ■ 0NU*C0NS2 

COMPUT39 

EOMU  > EOMU«CONS2 

COMPUTAO 

CE  • CE«CO*RMO 

COMPUTAl 

CE2  ■ C£2*C0*RM0 

COMPUTA2 

E20MO  ■'E2r)MU«C0NS3 

C0MPUTA3 

E30MU  ■ E3DMU«CONS3 

COMPUTAA 

ZWEMA  ■ ZWCMA«CUNSA 

computas 

Z2DMA  ■ Z2mMA*CON$A 

COMPUTAO 

COMPUTA7 

T1  • OA*X(3)»(TCRMl-TCRM2) 

C0MPUTA8 

TI  ■ T1  • ZZpMA  - EMAS 

C0MPUTA9 

T2  > EOMU 

COMPUTSO 

T3  • CE2 

C0MPUT51 

TA  ■ MA(NUM)«E(NUM) (UVNUM  • EZMAZ  * E30MU  - 

Z2WMA  * BMM 

C0MPUTS2 

NL  ■ Tl  ♦ T2  ♦ T3  ♦ TA  ♦ BMM 

C0MPUTS3 

T5  ■ MASS(X(3)*(TERM2-TEPM1) 

COMPUTSA 

T5  ■ T5  ♦ ZWMA  - ZMA 

COMPUT55 

TO  • -OMU 

COMPUTSO 

T7  ■ -CE 

COMPUT57 

111  CONTINUE 
RETURN 
ENO 

SUBROUTINE  COMPUT (X) 
dimension  X(6I 
REAL  RARtKoI 

REAL  NLtMASS«NCG«MtIT«IA*K,MAiMMAX«N 
INTEGER  END 
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BEST'AVAIUBIE  COPY 


c 

c 


c« 


TB  • •M«(Ni|M»*UVNUM  - eiOMU  * zmcha 
Bf  ■ BF/CXA 

fi.»TS*T6*TT.TB-Br 

ir(NPmNT,LT,3)00  TO  30 
IS  CONTINUE 

«RlTe<B»tO)NLtrU 

19  fORMATl"  NL  ■ "tEII,*,"  ft  ■ ••tCll*«) 

RETURN 

INO 

SUIROUTINC  INPUT 

• •••••  oeriNiTiON  or  input  vARUBtct 

XA  • INITIAC  TIME 

xE  ■ riNAi  tine 

HMIN  ■ NINIMUM  BTEP  SIZE 
HMAX  ■ NAX1NUM  STEP  SIZE 

EPSE  > RELATIVE  ERROR  CRITErIUN  USEO  POR  VALUES  OP  Y OT  A 
EPS  ■ ERROR  CRITERION  IN  RUTMER 

A ■ ABSOLUTE  ERROR  CRITERIA  USEO  IN  KUTNER 
NPrINT  • 1 PINAL  PRINTOUT 

■ 2 MATRIX  INVERSE  NATRlXfP  COLUMN  MATRlX»ANO  KUTMER 
RESULTS 

> 3 INTEORAL  VALUES 

> A CALCULATED  VALuES-CONSTANT  POR  OlVEN  INPUT  VALUES 
NPLOT  “0  NO  PLOT 

• I PRINTER  PLOT 
END  ■ NjMBEr  op  runs 


M ■ mass  op  CRAPT 
U ■ VEIiSHT  OP  CRAPT 
TZ  > THRUST  COMPONENT  IN  I UIrECTION 
Tx  • thrust  component  in  X OIrICTION 
xico  • distance  Prom  eo  to  center  op  pressure  por  normal  porce 

XP  ■ MOMPNT  ARM  OP  PROPELLER  THRUST 

XO  • DISTANCE  PROM  CC  tO  CENTER  OP  PRESSURE  POR  ORAO  PORCE 
kaiii*  aooeo  mass  COEPPICIENT 

AN  array  OIVEN  THE  VALUE  MAR  mHICH  IS  READ  IN 
BM(II  a BEAM  AT  PREE  SURPACE  UR  AT  CHINE 
ORAO  a PRICTION  ORAO 
K a MAVP  NUMBER 
RO  a HAVE  HEIGHT 
NU  a hAVP  StOPE 
NUM  a number  OP  STATIONS 
BL  a boat  LENOTH 
LAMBDA  a WAVE  LENOTH 

RO  a RADIUS  OP  OENERATION  IN  PEET 

T a Propelled  thrust  in  lbs 

OAHHA  a PROPELLER  THRUST  ANOLE  IN  UEOhEES 
OELTASaSTATlON  SPACINO  IN  PEET 

ECO  a LONOITUOINAL  CENTER  OP  ORAVlTV 
NCO  a vE-'TICAL  CO 
BETA) I)  a oEAO  rise 
NO(I)  a HEIOHT  OP  MEAN  BUTTOCK 

RHO  a density  op  hater 

ORAVTV  a OoAVlTY  PT/SEC*aZ 
OPR  <9  OCOoEES  PER  RADIAN 
RPD  a radians  per  OEOREE 
PI  a 3.1A159  


COMPUTSB 
C0HPUTS9 
COHPUT60 
C0MPUT61 
COHPUTBZ 
COMPUT63 
COMPUTBP 
C0MPUT6S 
C0MPUT66 
COMPUTOT 
C0HPUT6B 
INPUT  Z 
INPUT  3 
INPUT  A 
INPUT  ? 
INPUT  « 
INPUT  T 
INPUT  B 
INPUT  9 
INPUT  10 
INPUT  11 
INPUT  12 
INPUT  13 
INPUT  lA 
Input  is 
INPUT  16 
INPUT  17 
INPUT  IB 
INPUT  19 
INPUT  20 
INPUT  21 
INPUT  22 
iNPur  23 
INPUT  2A 
INPUT  25 
INPUT  26 
INPUT  27 
INPUT  2B 
INPUT  29 
INPUT  30 
INPUT  31 
INPUT  32 
Input  33 
INPUT  3A 
Input  35 
INPUT  35 
INPUT  37 
INPUT  38 
INPUT  39 
INPUT  AO 
INPUT  Al 
INPUT  A2 
Input  a3 

INPUT  AA 

input  as 
Input  a6 
Input  at 
Input  ab 

INPUT  A9 


oooooooooo 


BBT  AVAIUBIE  COPY 


c 

c 


c 


c 


c 

c 

c 


c 

c 


IST<1>  ■ STATION  POSITION 

iTAPT  ■ start  TIMC  or  THE  RAMP  ruNCTlON  FOR  SEA  wAVC 

RISE  ■ OURATION  or  THE  RISE  rHOM  ZERO  TO  ONE  Or  THE  NAHP 

••••tit  1C  OPTIONS 

IC(1)  ■!  OSE  NAVE  Z DISTANCE  IN  COHPUTINO  LlPT  COMPONENT 

or  NC  AND  ri 


REAL  lT«K«LAMBOA«M|MAiMMAX«NUfNtNCO«NOtMASS«NtilA,KAR 
INTEOCR  ENn 


COMMON  /const/  NC0«CCa«Pl«0PR»nPUt&MA)/TT«RH0»KiNUM«MA(l20)«CD»TAi 

• R(l20>fBETA«HH(iaO>«Ti«ORAO»R«XD»T«XPiM«IT« 

4 OCLTJkS<TX«EST(120)«CiRO«NAR»HMAX()''0)iTEST(120)* 

• N<lEO)iPHAir 

COMMON  /SHIP/  MASS«CINTiOA«CE«CE2fCE3«OMU«EOMOtE20MUiC)DMO»Br«BMM, 

• NL»rL«IA,C<120> 

COMMON  /IN/  BM(120l«ei(ll0>«vEt.|N 

COMMON  /INP/  N0(120) tXAiXEiHMAX»HMlN«A(6) »EPSE(b) iLAMBDA 

common/out/nprint«nplot«eno 

COMMON  /ACrlt/  XACCLtSwACCtCOACLiUL 

NAMCLIST/HSP/A«NPRINT«NPLOTiCN0,H«HL«TZ.TX4XCCa»XP«XD» 

• ORAOfROiTtOAMMAtECutNCOtKAR.ROttAMBDAiNUMtBETAiEST 

I tXAtXC»HMtN«HMAX,EPStVCLlN 

DATA  A /.01«,000l«.00001t.l«.u00l»»0000i/ 

DATA  NPR1NT«N^LOT<END/I«1iI/ 

DATA  *(,BI.,TZ*TX»XECO,XP»XO*OflAO, HO, lambda, RO,T»OAMMA, 

• ECOiNCOiKAR  /l6«,3«TS,6*0.0«*0A16,22,S,.9Sb2t2*0,0» 

t 2.32Si0*0il.0> 

DATA  NON, BETA, EST  /77,20.0, 

• 0,0000 «.0312S,,062S0«.ftO37S,.12S0O,.lS62S,. 16780, «21B7S, 

• ,2S060,,2012S,.312S0i.i<»S7S,. 37560,.  A062S,.07S6«.A667S, 

• .80000,. S3l2S,.802SO,.S037S,.b2S60,.bS62S,. 6875,. 71075, 

4 . 7S060,. 7012$,. S12S0,,8A37S,.B7S60,,9062S,, 93750,, 96875,1. 000, 

» 1.06250. I, 12500, 1.107*0, 1.25000, 1,3125, 1,37500, 1.9375, 

• 1. 500, 1. 5625, 1. 625, 1. 6675, 1.75il. 0125, I.b7S, 1.9375, 2.0, 

« 2. 0625, >.125, 2. 1975, 2. 25, 2. 3125, 2. 375, 2. 9375. 2. 5, 2, 5625, 2, 625, 

• 2. 6075, P.75. 2. 0125 ,2. 0750, 2. 9375,3.0,3. 0625, 3. 125, 3. 1975, 

• 3.2500,3.3125  ,3.375,3.9375,3.5,3.5625,3,625,3.6875,3.75  / 

DATA  XA,XC,MMlN,MMAX,r.PS  /O ,0,20,0, ,025, . 1 , , 1 5/ 

DATA  VELIN  /19,62/ 

••••••,  READ  IN  AND  WRITE  OuT  HuTMER  PARAMETERS  AND  PROGRAM 

OPTIONS 

REAOlS,NSP| 

WR1TE<6,HSP| 

DO  10  1*1,9 
10  EPSCU)  • *PS 

•••••••  SET  UP  constants 

PI  ■ 3.19I5926S3S99 
ORAVTY«32.lB 
OPR*57.2957T951308 
RP0«,0I7a532925I9 


NPUT  SO 
NPUT  51 
NPUT  52 
NPUT  33 
NPUT  59 
NPUT  SB 
NPUT  56 
NPUT  57 
NPUT  SB 
NPUT  50 
NPUT  *0 
NPUT  il 
NPUT  62 
NPUT  63 
NPUT  69 
NPUT  65 
NPUT  66 
NPUT  67 
NPUT  68 
NPUT  69 
NPUT  TO 
NPUT  71 
NPUT  72 
NPUT  73 
NPUT  79 
NPUT  75 
:nput  76 
:nPUT  77 
NPUT  78 
NPUT  79 

:nput  bo 

INPUT  SI 
INPUT  82 
INPUT  83 
INPUT  89 
INPUT  85 
INPUT  86 
INPUT  97 
INPUT  88 
INPUT  09 
INPUT  90 
INPUT  91 
Input  92 
INPUT  93 
INPUT  99 
INPUT  95 
INPUT  96 
INPUT  97 

Input  98 

INPUT  99 
INPUT  100 
INPUTlOl 
1NPUT102 
INPUT103 
I NPUT 109 
INPUTIOS 
INPUT106 
INPUTIOT 
INPUTIOB 


4 


1 


.i 


h 


05 


BESTAVAIHBIE'COPY 


If (EST(NUM|,CT.3,75»  STOP  3 
C 

C COMPUTE  NO  AND  AM  AARAVS 

C 

DO  32  l*liNUM 

irtesTiD.ac.o.TS)  oo  to  jo 

NO<l)>-0.4ASTi*(1.0*SaRT(EST(n/0.37S*<EST(n/0*75>*«2.0)) 
eM(I)a.37S«SORT(|.0*<EST(l)/.75-l*)««2.0I 
OO  TO  32 
39  NO(ll>0.0 
OMin  ■ 0,37S 
32  continue 

C«***»«*C0MRuTe  constants  ano  initialize  arrays 

MbM/ORAVTV 

RHOal.S'A 

IT«M*R0*R0 

K > 2 /LAMBDA 

CbSORT(ORAVTV/K) 

NUaRO*K 

PHALf  K <PI/2.)*RH0 
C 

BETA  ■ BETA*RPD 
CO  ■ COSIBfTAt 
TA  ■ TANEBETA) 

OO  60  laltNUM 

c<i>  > eco-esT(i) 

Nd)  ■ Nca*N0<n 
MMAl(I)  ■ KAR«PHALf*BM<n*BM(U 
TCST(I)  • (2.*BM(I)«TA>/PI 
69  continue 
ENOaCNO*! 

RETURN 

ENO 

SUBROUTINE  PLUTER(fX«XAiHMAX«LAMaUA*IBfNWAVE«lPT) 

C 


c 

INPUT* 

c 

Px 

A TmO  dimensional  ARnAY  containing  pitch  AND 

c 

HEAVE  VALUES  AT  EACH  TIME  STEP 

c 

XA 

INITIAL  TIME 

c 

MMAX 

TIME  INTERVAL,  PTIME*HMAX  • INTERVAL  BETWEEN 

c 

fX  VALUES 

c 

LAMBUA 

HAVELENOTH  USED  IN  CALCULATING  PITCH  AND 

c 

HEAVE  RATIOES 

c 

IS 

Number  or  rx  values 

c 

c 

nmavE 

START  Of  VALUES  AfTER  NAVE  IS  COMPLETELY  ON 

c 

REAL  iTfK 

,LAMaOA,M,MA,MMAX(N,NCO 

INTEGER  ENO 

c 

DIMENSION 

PX<2f600) ,fMIN(2>  tfMAX(2>,NVAR(2) 

C 

COMMON  /CONST/  NCOiCC9»PIiOPRtRPO*GflAVTy*RHOiK,NUM»MA(120l tCD*TAt 
k B(I2A| (BETAtHPdZOl tTZtORAO.M«XO(T*XP«MtITt 

• OELTASiTX,EST(120)«CtROtMA«MMAX<l20) iTESTdZOi* 

• Nd20|«PHALf 
COMMON/OUT /NPfl I NT  fNPLOT lENO 

C 

C •••••••  • SET  UP  VALUES  fOW  PLOT  AND  CREATE  PLOT 


1NPUTI09 
INPUTIIO 
iNPUTlIt 
INPUTllZ 
1NPUTII3 
INPUTIU 
INPUTIIS 
INPUTIIO 
INPUTII7 
INPUTIIB 
INPUTIIO 
INPUT 1 20 
INPUT121 
INPUT 122 
INPUT I 23 
INPUT  126 
INPUT  1 2S 
INPUT  126 
INPUT  127 
INPUT  128 
INPUT  129 
INPUT 130 
INPUT  131 
INPUT 1 32 
INPUT  133 
INPUT  136 
INPUT  135 
INPUT  136 
INPUT  1 37 
INPUT  136 
INPUT  139 
INPUT  160 
INPUT  161 
MLOTIR  2 
PLOTER  3 
PLOTER  6 
PLOTER  S 
PLOTER  6 
PLOTER  7 
PLOTER  6 
PLOTER  9 
PLOTER 10 
PLOTER 11 
PLOTER 12 
PLOTER13 
PL0TER16 
PLOTERIS 
PLOTER 16 
PLOTER17 
PLOTERIS 
PLOTER19 
PLOTER20 
PL0TER21 
PL0TER22 
PL0TER23 
PLOTER26 
PL0TER25 
HL0TER26 
PLOTER27 
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I 

i 


f 

r. 

t 


I. 

f; 


t 


f 


NfUH«2 

C ••»•••«•  SIT  UP  MIN  AKU  MAX  LIMITS  POM  PLOT 
rMtN(i)"rx(i»l> 
fMlN<2)«rX(2in 
FNAX(l»arX(ia) 

C • • UP  NIN  AKO  MAA  tIHIMTS  FOR  -ITCH  AND  MCAWC  RATIO 

FMNPaPXIZtNMAVl) 
rHAP>rx(Z*r4MAVC) 

FMNHaPxdtNWAVC) 

PMXHaPXdtNMAVO 


00  200  lal.IB  . , . 

ir(fxU«i).LT.rHiHdnrMiNd»-Fxj  ti) 

ir(Fxd#n.0T.rMAxddfMAxmaF,d,i) 
IF<FX<t*n.LT.F«lN<2))FMlN(2>afK(Z,n 
lF(FX(2»n.0T.FMAXI2»|FMAX<|»arx(2*n 
IFd*Le«HHAVC>00  TO  200 
IF (FX d ♦ I ) .UT .FMNHJFRNMaFX d , 1 » 
IF(FXd*n.OT.FMXH»FHXMaFXd»n 
IF <FX(2* I >.LT.FMNR)FMNRarX(2fl) 
IF<FX<2d».0T.FMXR»FMXRaFX(2»I» 

ZOO  CONTINUE 

IFdPT*€O«0»  00  TO  600 
C»»#aaaa  COMPOTC  RAT  IOCS 
cots  a (FMXH-FHNH)/<2**N0) 

COLA  a <rMXP-FHNP>/dA.aPI»RO»/LAMiOA> 
RRITC(A»T00)  C0L3»C0LA 
TOO  FORMAT (IMI,"  MCAVC  AMPLIIUOC/RAVfcMtlOMT 


•SC12tO*/«ZXi 


PdCM  lMPLlTU0l/«2.iPI*WAVtHtl0hT/LAMB0AI  a ...CIZ.O) 


000  CONTINUE 

■ NVARd)*IOH  MCAVI 
NVARI2»*I6m  FITCM 

N|a2 

X0«XA 

IF'tNp’oT^Eg.nCALL  PLOT2tFX,rNIN|FMAX,NVAR,NFuN,Nl,Ii,XO«OeLX> 

rcturn 

ENO 

SUBROUTINE  TRAP(F*DX»NPTStANS) 


input  I 

F 

OX 

NPTS 

UUTPUTI 

ANS 


ARRAY  OF  FUNCTIONAL  VALUES  OF  THE  INTEORANO 
THE  X INTERVAL  BETWEEN  VALUES 
THE  NUMBER  OF  VALUES  OlVEN 

THE  value  of  THE  INTEGRAL 


DIMENSION  F(NPTS) 

ANfaO.O 

IF<NPTS»LT.Z)00  TO  900 
00  I lalfNPTS 

I ANSaANS*Fd)  , , 

ANSaOXa(AN5-0.5*<Fd»«F(NPT9»  I > 


999  CONTINUE 
RETURN 

END  . . 

FUNCTION  RMPdtSTARTtRlSE) 


HLOTERZe 

PL0TERZ9 

HLOTEP36 

ML0TEP3I 

PL0TER3Z 

PLOTER33 

PL0TER3A 

PL0TER3S 

PL0TER36 

PL0TEB3T 

PLOTCB30 

PL0TER39 

PLOTERAO 

PLOTERAl 

PL0TERA2 

PL0TCRA3 

plotcbaa 

PLOTCRAS 

PLOTERAO 

PL0TERA7 

PLOTERAO 

PLOTERAO 

ploterbo 

PLOTERSl 

PL0TERS2 

PL0TERS3 

PL0TER8A 

PLOTERSS 

PL0TERS6 

PLOTERST 

PLOTERSS 

PLOTERBO 

PLOTER60 

MLOTEROI 

PLOTER02 

PL0TER03 

PLOTEROA 

PLOTEROS 

PLOTEROA 

PLOTEROT 


TRAP 

2 

TRAP 

3 

TRAP 

A 

TRAP 

3 

TRAP 

0 

TRAP 

7 

TRAP 

0 

TRAP 

9 

TRAP 

10 

TRAP 

d 

TRAP 

IZ 

TRAP 

13 

TRAP 

lA 

TRAP 

19 

TRAP 

10 

TRAP 

17 

TRAP 

18 

TRAP 

19 

RMP 

Z 
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This  FUNCTION  IS 

USED  TU 

ChADUALLY  IHPLIHENT  ThE  WAVE 

i<MP 

3 

kMP 

A 

i 

T 

CURRENT  TIME 

MNP 

s 

•t 

START 

TINE  TO  START 

ramp  pru»< 

1 0.0 

TO  i.O 

MMP 

6 

/ 

Rise 

THE  LENOTM  OR 

THE  rue 

fROM 

0.0  TO  1.0 

RMP 

7 

PMP 

0 

1 

H*0»0 

HMP 

9 

1 

ir<T.LT.START>00  TO  99 

RMP 

10 

1 

lf(RlSBteO. 

0.0)00  TO  SO 

RMP 

11 

' 

top«t-start 

RMP 

12 

H«1*0 

RMP 

13 

ir <tor.lt. RiieiM>TOR/Rise 

RMP 

19 

C 

00  TO  99 

HMP 

IS 

! ^ 

•0 

H>1. 

RMP 

10 

ircT.eo.STARnn-o.s 

RMP 

17 

♦9 

RMR«M 

RMP 

10 

; 

RETURN 

RMP 

19 

CNO 

HMP 

20 

I 


i 


;3 


3 

•1 
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LISTING  OF  COMPUTER  PROGRAM  FOR  CALCOMP  PLOTS 


PROGRAM  PLTHSP (INPUT » OUTPUT . TAPES- 1 NPUT t T APES-OUTPUT » T APE T , T APE« > 
ITAPE  « r 

CALL  CALPLTdTArCl 

STOP 

END 

SUBROUTINE  CALPLT(tTAPC) 

UIMENSIUN  time (A003) tPITCHiABOSt  *nEAVE (40031 
• tIBUPdOOO)  tBHACL  (4043)  tCGACL  (4003) 

LOGICAL  ACCEL 

CAL  CUMP  PLOT  OP  PITCH  ANO  HEAVE  VERSUS  TIME 
IREAO  ■ S 

REAOdREAOtlO)  XAAlS«rAXlSP*rAAtSN«MT 
10  FORMAT  ( (IP  10.0) 

' ACCEL  > tP^LSE. 

REAOdREAO.SO)  lA 
20  PORMATdlO) 

' IPdA.EO.l)  ACCEL  ■ .TRUE. 

IPIACCELI  oEAUilREAOtlOl  VAxlSt),yAAlsC 
CALL  REAUT (TIME (HEAVE tPITCH.BDACLtCGACLtNPTbl 
PLUTS(iauP«1000«7) 

PLUT(0,8*l.9t-3) 

ESCALC(TlMe.XAXlS«NPTS«l> 

ESCALP (HEAVE • VAX 1SHiSPT8< 1 ) 

ESCALP(PlTCH|YAXlSP«NPT8d) 

CALL  CSCALC<BHACL*YAxlSH*NPTSd) 

CALL  ESCALEICOACLtVAXlSC.NPTSd) 


CALL 
CALL 
CALL 
CALL 
CALL 

IP (ACCEL) 

IP (ACCEL) 

N1  • NPTS*) 

N2  ■ NPTS*P 
N3  • NPTS*3 

CALL  CAXlS(O,O,0.OdSHTlME  IN  SECUNDSi-lS.XAXlS.O.O* 

« TIMelNl)dlMC(N2)»TlME(N3)»HT) 

CALL  CAX|S(0«0«O.Od3HHEAVE  IN  PEtTd3tYAXISH*V0*0« 

* HEAVE (N|) ;HeAVC(N2) f HEAVE <N3> 'HT) 

TEMP  ■ TIME(N2) 

T1ME(N2)  ■ T1ME(N2)/T1HC(N3) 

HEAVE (N2)  ■ HEAVE (N2) /HEAVE (N3) 

CALL  LINE(TIME«HEAVE»NPTStl<OtO) 

TIMC(N|)  ■ TEMP 
XNEW  ■ XAXISO. 

YNEN  ■ 1.0 

CALL  PLUT(XNC»(iO.Of3) 

CALL  EAXIS(O.O,0.Od5HTlME  IN  SECUNUS.-lS.XAXlS.O.O* 

* TlMC(Nlj .TIME(N2).T1ME(N3),mT) 

CAI.L  EAXlS(O.O,0.Od3HPlTCH  IN  RAU.dS. YAXlSPfVO.O. 

, P1TCH(NI  j .PITCM{N2)  .PlTCn(N.T(  tHT) 

TIME(N2)  ■ TIME(N2)/TIMC(N3) 

PITCM(N2)  • P1TCH(N2)/PITCM(N3( 

CALL  LINC(TIHt,PITCrt.NPTS(l»0»0) 

IP (.NUT. ACCEL)  GO  TO  30 
TIMC(N2)  ■ TEMP 
CALL  PLUT(XNC«i0.0.>3) 

CALL  iAXlSIO.OtO.O.lSHTlME  IN  StLUNOS.-IS.XAXlSf O.O.T IME (Nl ) « 

. T1hE(N?) ,TIMC(N3) iHT) 

CALL  LAXIS0.0.0.0.16HRUW  ACCELEHAT iUN. 1 6, YAXl SB.90 .0 iHHACL (Nl ) 

* dwACL(N2)  fB(|ACL(Nj)tHT) 

T1ME(N2)  ■ TlMC(Ne)/TlHE(N3) 
dVACL(N2t  - UnACL(N2)/BmACL(N3) 


2 

3 

4 

5 

6 
2 

3 

4 

s 

6 

r 

a 

9 


main 
MAIN 
MAIN 

main 

MAIN 
CALP 
CALP 
CALP 
CALP 
CALP 
CALP 
CALP 
CALP 
CALP  10 
CALP  1 1 
CALP  12 
CALP  13 
CALP  14 
CALP  IS 
CALP  16 
CALP  17 
CALP  IB 
CALP  19 
CALP  20 
CALP  21 
CALP  22 
CALP  23 
CALP  24 
CALP  2S 
CALP  26 
CALP  27 
CALP  20 
CALP  29 
CALP  30 
CALP  31 
CALP  32 
CALP  33 
CALP  34 
CALP  3S 
CALP  36 
CALP  37 
CALP  30 
CALP  39 
CALP  40 
CALP  41 
CALP  42 
CALP  43 
CALP  4<. 
CALP  45 
CALP  46 
CALP  47 
CALP  48 
CALP  49 
CALP  50 
CALP  51 
CALP  52 
CALP  53 
CALP  54 
CALP  55 
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CALL  LINe(TIML«t)WACL*NPTSf  1>0,0» 

TlMe^N^^  ■ TCMP 
CALL  PLUT(XNeHi0*0»-3l 

CALL  CAXtSIO.OtO.OilShTlMC  IN  SfcCUNDSt-lStXAXlStO.OtTlML INI) • 

• TIPeiNS) iTIMElNS) thT) 

CALL  LAXIS(0.O,O,0ilSHCl'  ACCeLex(ATII>N(lStVAxlSCt90.0«COACL (N1 ) • 

• C0ACL(N2)fC0ACLIN3)tHT) 

T1ME(N2)  ■ TIHC<N2)/TIHE<N3) 

C0ACL<N2)  ■ C(tACL(N2)/C0ACL(N3) 

CALL  LlNCalMCtCaACLtNPTStltOtO) 

30  continue 

CAUL  PLOT  no. OiO. 0.999) 

RETURN 

END 

SURROUTINE  react (TlME.HCAvE.RITChtBMACL.COACL.NPTS) 

DIMENSION  X(6)  .HEAVEa)  iPITCHd) 

• .TIMEd)  iRMACL(I)  iCOACLd) 

1 • 0 

S CONTINUE 
1 > 1*1 

READ (9)  TlMEd)  «(Xd)«I*A.6),UNACLd).CbACL41) 

1P(EOPI9))10.1S 
15  CONTINUE 

WRlTC<6i20l  TIMEd) t (X(JI .JbA.A) .dWACLt I) .COACL (I ) 

29  FORHATdH  .6(F7.2.2X)) 

HEAVEd)  ■ X(S) 

PlTCHd)  ■ X(«) 

Ifd.OE.AOnO)  00  TO  10 
00  TO  5 
10  CONTINUE 
■ NPTS  • 1-1 
RETURN 
END 

SUBROUTINE  EAXlSlXPAOE.YPAOEf iBCO.NCHAR.AXLEN.ANOLC.riRSTV. 

4 OEUTAV.OELTAU.HT) 

DIMENSION  tBCOd) 

THIS  RUMTlNE  WORKS  LIKE  THt  CALCUMP  AXIS  WITH  THE 
EXCEPTION  THAT  THE  TICK  HANKS  AnE  NOT  NECCESSARILY 
EVERY  INCH  AND  THE  HEIOht  OF  T^E  CHARACTERS  IS  INPUTTEO 

CAUL  PLOT(XPA0E.YPAaEi3) 

ISN  ■ ISIONd.NCHAR) 

ISON  ■ SiaNd..OeLTAV) 

AMIN  ■ riRSTV 
X > XPACE 
V ■ rPAOE 
XNUH  ■ fIrstv-oeltav 
N ■ AXLEN/OELTAU 
IF(N*0FLTAU.LT.AXLENI  N-N-l 
AHAX  a AMIN*(N«0ELTAV) 

NOlO  a NDlOITIAHIN.AHAx.OLLTAUtNO) 

10  CONTINUE 

TEST  a (NOIOaHT)  ♦ MT 
IFdESY.OT.UELTAU)  HTahT/2. 

IFITFbT.OT.OELTAU)  00  TO  |0 
AYN  a d.SaHT) 

BYN  a ( ( (N0ia-2)aHT)/2.*.bahT) 


CALP 

S6 

CALP 

57 

CALP 

SB 

CALP 

59 

CALP 

60 

CALP 

61 

CALP 

62 

CALP 

63 

CALP 

6A 

CALP 

6S 

CALP 

66 

CALP 

67 

CALP 

68 

CALP 

69 

CALP 

70 

HEAD 

2 

HEAD 

3 

READ 

4 

READ 

s 

READ 

6 

READ 

7 

READ 

B 

READ 

9 

READ 

10 

READ 

U 

READ 

12 

READ 

13 

READ 

U 

READ 

15 

READ 

16 

READ 

17 

READ 

IB 

READ 

19 

READ 

20 

eaxIs 

2 

eaxis 

3 

EAXIS 

A 

EAXIS 

5 

EAXlS 

6 

EAXlS 

7 

EAXIS 

S 

EAXIS 

9 

EAXlS 

10 

EAXIS 

11 

EAXlS 

12 

EAXIS 

13 

EAXIS 

14 

EAXlS 

15 

EAXlS 

16 

EAXlS 

17 

EAXlS 

18 

EAXIS 

19 

EAXIS 

20 

EAXIS 

21 

EAXlS 

22 

EAXIS 

23 

EAXlS 

24 

EAXlS 

25 

EAXIS 

26 
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N • N*l 

LAXIS 

27 

TANG  * (QO. •anglo/57. 2958 

bAXiS 

28 

ANQ  > ANGLt/57.?954 

bAKl3 

29 

5T  ■ SINITANO) 

eaxIs 

30 

CT  • COSITANO) 

eaxis 

31 

S ■ STN(ANU) 

EAXIS 

32 

C ■ C0S<ANG) 

EAXIS 

33 

00  30  1«1«N 

Eaxis 

3A 

ird.LO.l)  bU  TO  20 

EAXiS 

3S 

A ■ k»oi:ltau*c 

eaxis 

36 

t ■ y»ocltau*$ 

EAXlS 

37 

CALL  >'LUT(X*t«2l 

EAXIS 

30 

ir(UCO.N)  OU  TO  20 

EAXIS 

39 

XT  • X*(,1*CT«I»NI 

EAXIS 

40 

YT  ■ YM.I*ST*ISNi 

eaxis 

41 

CALL  ►*L0T(XT#YT.?) 

EAXiS 

42 

2^ 

XN  ■ X«AYN*CT*1SN-AYN«C 

EAXlS 

43 

YN  ■ Y*AYN»5T»I5N-BYN*S 

eaxis 

44 

XNUH  ■ XNUM*D1LTAV 

eaxis 

4S 

CALL  NUMe(;RtXN«YNtHT«XNUMtANGLL*NO) 

EAXIS 

46 

CALL  ►*L0T(X»Yf3) 

eaxis 

47 

39 

CONTINUL 

EAXIS 

40 

XAB  ■ ( ( (AxLEN/HTt/2.)«(IABS(NCHAH)/2.M*HT 

Eaxis 

49 

YSO  • 3,8«hT 

eaxis 

SO 

XT  • XPAGt  ♦ XS*»*C  ♦ t$N«YSf'*CT 

Eaxis 

SI 

YT  ■ Yt»Aat  • XSP«S  ♦ isn^ysh^st 

Eaxis 

92 

CALL  SYMBOL ( X T « YT  «H  r » 1AC0« ANGLE • 1 AdS ( NCHAA ) ) 

Eaxis 

53 

RCTURN 

EAXIS 

94 

ENO 

EAXlS 

SS 

Function  NoioiT(AHiNtAMAx»ANuN»NU) 

NDIO 

2 

NOlO 

3 

riNos  thc  number  or  digits  necclssary 

TO 

PRINT 

NDIO 

4 

EVEN  increment  or  THE  rUNCTlON  ON  THE 

AXIS 

NDIO 

9 

NDIO 

6 

NOIOIT  THE  NUMBER  or  PLACES  IN  THE 

ENTIRE 

NUMBER 

NDIO 

7 

NO  THE  NUMBER  Or  DECIMAL  PLACES 

NDIO 

0 

ANUM  THE  VALUE  OtVEN  TO  EACH  INCREMENT 

ON  THE  AXIS 

NDIO 

9 

NDIO 

10 

ir(ABS(AMlN|,LT.ABS(AMAX))  00  TO  20 

NDIO 

11 

ir(AMS(AMlN),ea.ABS(AMAXUANO*AMAX.Ne.O) 

00 

TO 

20 

NDIO 

12 

ir(A0S(AMIN|.GT.ABS(AMAX>)  00  TO  10 

NDIO 

13 

AMAX  ■ 1, 

NOIC- 

14 

AMIN  • -I. 

NOIO 

IS 

00  TO  20 

NDIO 

16 

10 

AMAX  ■ ABS(AMIN) 

NDIO 

17 

20 

IFIAMAX.LC.l.)  00  TO  SO 

NDIO 

10 

NOIV  ■ in 

NDiO 

19 

I • 1 

NOIO 

20 

30 

IF<AMA*/N0IV,LT,1»  00  TO  ‘•0 

NDIO 

21 

I ■ I*I 

NOIO 

22 

NOIv  • NOIV*lO 

NOIO 

23 

00  To  no 

NOIO 

24 

40 

NOIOIT  ■ [*3 

NOIO 

25 

NO  a 2 

NOIO 

26 

00  TO  8U 

NDIO 

27 

S9 

NOIv  a lO 

NOIO 

20 

1 a 1 

NOIO 

29 

OO 

iriAMAxaNOlV.OT.l.l  00  TO  70 

NDIO 

30 

I a 1*1 

NDIO 

31 

71 
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NDIV  • Nt)IV*10 

NOIO 

32 

00  TO  60 

NDIO 

33 

TO  NOIOIT  ■ 1*2 

NOIO 

3A 

' NO  ■ 1 

NOIO 

3S 

ao  00  ■ rLUAKNOt 

NOIO 

36 

A ■ ANUH«(10»*00) 

NOIO 

3T 

IX  ■ X 

NOIO 

3a 

ir(X-FLUAT(IX>.LT..0001)  00  TU  90 

NOIO 

39 

00  ■ 00*1 

NOIO 

AO 

NO  ■ N0*l 

NOIO 

A1 

NOIOIT  ■ NOIOIT*! 

NOIO 

A2 

00  TO  ao 

NOIO 

A3 

90  CONTiNut 

NOIO 

AA 

aCTUAN 

NOIO 

AS 

CNO 

NOIO 

A6 

SUaaOUTlNC  CSCALC  < APKAV .AXLCNtNATSi inc> 

ESCAL 

2 

escal 

3 

riNOs  thi  scale  to  at  used  on  toe  axis  • 

ESCAL 

A 

ARNAT  MiiST  HAS  THREE  UNUSEO  POSITIONS 

ESCAL 

S 

ARPAT(NPTS*n  ■ PIRSTV 

ESCAL 

6 

ARoAV(NPTS*2>  > OELTAV  <ThE  INCREMENT  OETPEEN  TICK 

MARKS  escal 

7 

VALUES  • NUMBERS) 

ESCAL 

a 

ARoAV<NPTS*3l  • OELTAU  (THE  INCREMENT  IN  INCHES 

ESCAL 

9 

METBEEN  TICK  HARKS  ) 

ESCAL 

10 

escal 

11 

ESCAL 

12 

OIMENSIUN  array <1) 

ESCAL 

13 

AMIN  ■ ARRAV(I) 

ESCAL 

lA 

AMAX  ■ array <1) 

escal 

15 

ISON  ■ iSIONtUlNC) 

ESCAL 

16 

INC  • iabsunci 

ESCAL 

17 

OU  10  Ial«NPTS«INC 

ESCAL 

IS 

If ( AORAY ( 1 ) .LT « AMIN)  AMIN-AHRAY 4 1 ) 

escal 

19 

If ( APRAY ( 1 > .OT.AMAX)  AMAXaAHHAY  4 1) 

ESCAL 

20 

1(1  CONTINitC 

escal 

21 

aO  AUNIT  ■ llNlT(AMtNtAMAXtAXLCNiNtANUM) 

ESCAL 

22 

CALL  AOUUAT<AMlN,AMAX«AuNlTtAXLCN|N,ANuM) 

ESCAL 

23 

ARRAY (NMrS«l)  ■ AMIN 

ESCAL 

2A 

ARRAY 4NHTS*2)  « ANUM*IS0N 

escal 

29 

ir(lS(iN,AO.-l)ARRAY4NPT8*l)  ■ AMAX 

escal 

26 

ARRAY(NHTS«3)  ■ AUNIT 

escal 

27 

If 4A8S4ANU'*) «C0. AUNIT)  ARRAY (NPTS*21  ■ 1,»IS0N 

ESCAL 

2S 

If (ARS(AMUm) .EO. AUNIT)  ARRAY (NP rS*3)  ■ 1. 

ESCAL 

29 

RETURN 

ESCAL 

30 

END 

escal 

31 

SUBROUTINE  ADJUST (AMIN« AHAX • AUNII tAALENiNtANUM) 

JUST 

2 

JUST 

3 

OIVEN  A'UN  ANO  AMAX  WHICH  ARC  UlSTINCT  VALUESt  ADJUST 

JUST 

A 

THEM  SO  that  they  are  even  MULTIPLES  Of  AUNIT 

JUST 

9 

JUST 

6 

K ■ 1 

JUST 

7 

MIN  ■ AMIN/ANUM 

JUST 

6 

If  4AMIN.lt, MIN»ANUM)  MIN  ■ MlN-1 

JUST 

9 

AMIN  ■ MIN«ANUM 

JUST 

10 

max  ■ AMAX/ANUM 

JUST 

11 

If 4AMAX.UT.MAX»ANUM)  MAX  ■ MAX*) 

JUST 

12 

AMAX  ■ MAX*ANUM 

JUST 

13 

10  TERM  ■ AMIN. (N«K)*ANUM 

JUST 

lA 

If (TERM. LT. AMAX)  00  TO  20 

JUST 

19 
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K ■ K«l 
00  TO  10 

20  AUNIT  ■'AAlCN/<N-K«I) 

N u AXLtN/AUNlT*l 
RCTliRN 
‘CNO 

Function  UNlT(AHlN,AMAX«AKLeN»N,ANuN> 

finds  the  INCMMCNT  BCTMCCN  values  TU  be  used  UN  THE 
AXIS  IN  AS  FAR  AS  LABELlNO  THE  TICK  MARKS 
FINOS  THE  NUMBER  OF  DIVISIONS  TO  HE  MADE  ON  THE  AXIS 
FINOS  THE  SHE  IN  INCHES  OF  THESE  DIVISIONS 

IF(AMIN.NE.AMAX)  00  TO  10 
AMIN  • AMIN-I 
AMAX  ■ AMAX*I 

10  lF(AMAX.LT,l.AN0,AMIN,01.«l)ao  TO  110 
30  MIN  ■ AMIN 
MAX  a AMAX 

IF (AMAX. OT. MAX)  MAXaMAX*! 

IF(AMIN.LT.HIN)  MINaMIN-1 
1F<MIN»LT.0)  NhIO  ■ MAX*tABS<MIN) 
lF(MIN,OC.n)  NhIO  a MAX-MIN 
NUN  a 1(1 

«U  IF(NHIO.LT«NUM)  00  TO  60 
NUM  a NUH*10 
00  TO  AO 

60  N a NHlO/TNUM/10) 

IF(MlNtCT40«ANO.MAX.OT«0)  00  TO  70 
IF(Na(NUN/10).UT.NHtO)  NaN*l 
ANUM  a NUM/10. 

AUNIT  a AXLEN/N 
00  TU  lAO 

70  NN  a IABS<HIN>/(NUM/10I 

IF(NN*(NUN/10>.LT.IABS<N1NM  NN  a NN*1 
N a MAX/<NUM/10) 

IF(N*<NUM/10) tLt.MAX)  N a N*1 

N a N«NN 

ANUM  a NUN/IO. 

AUNIT  a AXLEN/N 
00  TO  160 
110  NuMalO 

llO  IF(AMAxaNUH.OT.l)  00  TO  130 
NUM  a hum# IQ 
00  TO  120 

130  UNITT  a I, /NUM 

1a6  Nl  a AMlNaNUM 

N2  a XHAXaNUM 

IF(AMlNaNUM,LT,Nl)  NlaN|-l 
If (AMAxaNUNtOT.NS)  N2aN2*l 
IF(N1.NE.N2)  00  TO  ISO 
AMIN  a AMIN-UNITT 
AMAX  a AMAX-UNITT 
OU  TO  lAO 

ISO  N a NZ'Nl 
ANUM  a UNITT 

IF (AMIn.LT.O.ANO.AMAX.LT.O)  NaNl-NZ 
IF (AMlN«LT.OtANOtAMAX*OC«U)  NaNE-Nl 
AUNIT  ■ AXLEN/N 


OUST 

16 

JUST 

17 

JUST 

la 

JUST 

19 

JUST 

20 

JUST 

21 

UNIT 

2 

UNIT 

3 

UNIT 

h 

UNIT 

S 

UNIT 

6 

UNIT 

7 

UNIT 

8 

UNIT 

9 

UNIT 

10 

UNIT 

11 

UNIT 

12 

UNIT 

13 

UNIT 

lA 

UNIT 

15 

UNIT 

16 

UNIT 

17 

UNIT 

10 

UNIT 

19 

UNIT 

20 

UNIT 

21 

unit 

22 

unit 

23 

UNIT 

26 
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NICAL  VALUE  THEY  CARRY  A CONSECUTIVE  NUMERICAL  IDENTIFICATION  NfiOAROLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIQINATINO  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A PRELIM- 
INARY, TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
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BASIS. 


